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Abstract. We consider the 'classical' Boussinesq system in one space dimension and its symmetric analog. 
These systems model two-way propagation of nonlinear, dispersive long waves of small amplitude on the 
C " 3 ' surface of an ideal fluid in a uniform horizontal channel. We discretize an initial-boundary-value problem 

^N) ' for these systems in space using Galerkin-finite element methods and prove error estimates for the resulting 

semidiscrete problems and also for their fully discrete analogs effected by explicit Runge-Kutta time-stepping 
procedures. The theoretical orders of convergence obtained are consistent with the results of numerical 
experiments that are also presented. 
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1. Introduction 

In this paper we will analyze Galerkin approximations to the so-called 'classical' Boussinesq system 



< 

in 

(N 

< 

J2 ' which is an approximation of the two-dimensional Euler equations of water-wave theory that models two-way 

propagation of long waves of small amplitude on the surface of an ideal fluid in a uniform horizontal channel 
of finite depth. The variables in (ll.ip are nondimensional and unsealed; x and t are proportional to position 
along the channel and time, respectively, and 77 = 77(0;, t) and u = u{x, t) are proportional to the elevation 
^ ' of the free surface above a level of rest represented hy rj — 0, and to the depth-averaged mean horizontal 

00 I velocity of the fluid. 

\l ' The system (jl.ip is a member of a general family of Boussinesq systems derived in llj that are approxi- 

mations to the Euler equations of the same order as (jl.ip and whose nonlinear and dispersive terms are of 
equal importance when written in scaled form. These systems are written as 

(3 ■ 'nt + Uj; + (rjUjx + au^xx - or]xxt = , 

O . Mt + + UUx + cri^xx - duxxt = , 

where a, b, c, d are real parameters satisfying a + b — ^{9'^ — 1/3), c + d = i(l — 9^), where < 9 < I. 
. ^ , The specific system p.ip has been previously formally derived from the Euler equations, in the appropriate 

' parameter regime, in [53], [55], and |3T]. It has been widely used in the engineering fluid mechanics literature 

^ I for computations of long, nonlinear dispersive waves, starting with [24] and [25'. 

Existence and uniqueness of the initial- value problem for (jl.ip posed for x G R and < > and supplemented 
by initial conditions of the form 

r/{x,0) ~ r]o{x) , u{x,0) — uq{x) , x G K , (1-3) 

was studied by Schonbeck, ^7} and Amick, 2 . In these papers global existence and uniqueness was es- 
tablished for infinitely differentiable initial data of compact support such that 770 (x) > — 1, a; € R. In 
[l2] the theory of [27] and [2] was used to prove that given initial data {r]Q,uo) G H'^{R) x if'*+^(M) 
for s > 1, such that mfx^RijQ^x) > —1, there is a unique solution (7?,u) which, for any T > 0, lies in 
C(0,r;i?"(R)) X C(0,T;iJ^+i(M)). (Here, iJ^(R) is the usual, L^-based Sobolev space of functions on R 
and C(0, T; X) denotes the space of functions (p — 4'i't) having, for each t S [0, T], values in a Banach space 
X and are such that the map [0,r] 1— ||0||x is continuous.) 
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It is well known that the initial- value problem (I1.3P has classical solitary-wave solutions. In [3] 

we review the properties of such solutions, we construct them numerically and investigate, by means of 
numerical experiments, the resolution of general initial profiles into sequences of solitary waves, the details 
of the interaction of solitary waves during head-on and overtaking collisions, the interaction of solitary 
waves with boundaries, and several other issues. These numerical experiments were performed in |3] by 
fully discrete Galerkin methods that approximate solutions of p.ip posed on finite intervals and subject 
to boundary conditions. It is of interest therefore to study initial-boundary-value problems for and 
establish error estimates for their numerical approximations. 

In this paper we shall analyze the numerical solution of the following initial-boundary-value problem 
(ibvp) for (jl.ip : For some < T < oo, we seek rj — rj{x,t), u — u{x,t), defined for < .t < 1, < t < T, 
and satisfying 



This ibvp (for < t < oo) has been studied by Adamy, [Tj, who showed that it has weak (distributional) 
solutions {ijju) e _L°°(M+;L^ x Hq) provided e.g. that rjo S L^, uq € Hq with infa;g[o,i] r?o(a;) > —1. (Here 
= L^{Q, 1), and Hq = Hq{Q, 1) is the subspace of the Sobolev space H^{0, 1) whose elements vanish at 
X = and x — 1.) The proof uses a parabolic regularization of the first pde of (jCB[) . a technique used in 
the context of the Cauchy problem by Schonbeck, [27] . It should be noted that the homogeneous Dirichlet 
boundary conditions on u in (jCBp are one kind of boundary conditions that lead to well posed ibvp's in 
the case of the linearized system, [2l]. It is noteworthy that the CB system needs only two boundary 
conditions for well-posedness as opposed to the four boundary conditions (for example, Dirichlet conditions 
on 1] and u at each endpoint of the interval) required in the case of other Boussinesq systems, such as the 
BBM-BBM (a = c = 0, & = d in pT2]) .[T0]). or the Bona-Smith systems (a = 0, 6 = d > 0, c < in 
(ll.2p . 7 .) The case of the homogeneous Dirichlet boundary conditions in (jCB[) may be viewed as a base for 
studying the nonhomogeneous analog wherein u{0,t) and u{l,t) are given functions of t > corresponding 
to measurements of the velocity variable at two points along the channel. 

In [IB], Bona, Colin, and Lannes introduced another type of Boussinesq systems that they called 'com- 
pletely symmetric'. These are obtained by a nonlinear change of variables from the usual systems (|1.2I) . and 
have certain mathematical and modelling advantages over the latter. In this paper we shall also therefore 
consider the analogous problem for the symmetric system, specifically the ibvp 



for t > 0, which simplifies the study of its well-posedness and the analysis of its numerical approximations 
as will be seen in the sequel. 

In this paper we shall analyze semidiscrete and fully discrete Galerkin-finite element methods for the ibvp's 
(ICBI) and (jSCB[) assuming that their solutions are sufficiently smooth. Previously, rigorous error estimates 
for Galerkin methods for these 'classical' Boussinesq systems were proved in the case of the periodic initial- 
value problem in j3j and |9j . (For numerical work for this type and for other Boussinesq systems of the form 
(HID we refer the reader e.g. to [3], [5], [S], [22], M , [H], [13], [13, [H], [i ) 

The analysis of Galerkin methods for (ICBI) or (jSCBI) is of considerable interest due to the loss of optimal 
order of accuracy that emerges from the error estimates and is also supported by the numerical experiments. 
We shall investigate this phenomenon in detail in the paper but one may in general say that in the uniform 
spatial mesh case the (limited) loss of accuracy seems to be due to a combination of effects related to the 



Vt + u, + {m)x = , {x,t)e [0, 1] X [0, T] , 

Ut + T^x + uux - \uxxt = , (a;, t) e [0, 1] x [0, T] , 
77(2;, 0) = 7^0 (x) , u(a;, 0) = uo(a;) , a; G [0, 1] , 
w(0,t) = 0, u(l,t)=0, te[0,T]. 



(CB) 



m+Ux + \(tiu)x - , (x, t) G [0, 1] X [0, T] 

Ut + T^x + ^uux + \rjrjx - \uxxt = , {x,t) e [0, 1] x [0, T 
77(a;,0) = 77o(a;) , u{x,0) ^ uq{x) , a; e [0, 1] , 
M(0,t) = 0, u{l,t) = 0, t£[0,T]. 



(SCB) 




(1.4) 
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boundary conditions and the specific form of the 'classical' Boussinesq systems. (For example, no loss of 
accuracy occurs in the case of periodic boundary conditions on rj and u for these systems, cf. [3], [9].) In 
the case of general (quasiuniform) mesh the (more severe) loss of accuracy seems to stem from the lack of 
cancellation effects that are present in the uniform mesh case, and from the hyperbolic character of the first 
p.d.e. of these systems. (The loss of optimal order of accuracy in standard Galerkin semidiscretizations of 
first-order hyperbolic equations manifested in various contexts is well known and was early observed e.g. by 
Dupont in pD].) 

In section 2 we consider the standard Galerkin semidiscretizations of (|CBp and (|SCBI) in the space of 
piecewise linear, continuous functions. In the case of a quasiuniform mesh, we prove that if the solution (77, u) 
of (|CBp or (jSCB|) is sufficiently smooth and {rih,Uh) is the semidiscrete approximation, then ||77--?7;i|| = 0{h) 
and |ju — Uh\\i — 0{h). (Here || • ||, || • ||i denote, respectively, the and norms on [0, 1].) In the case 
of uniform mesh, a suitable superaccuracy result for the error of the interpolant into the finite element 
subspace (a consequence of cancellations due to the uniform mesh) affords proving the improved estimates 
Wv ~ VhW = 0{h^^^) and ||u — = 0{h?). These rates of convergence are confirmed by the numerical 
experiments at the end of the section. 

In section 3 we turn to the semidiscretization of (jCB[) or (jSCB[) using as finite element subspace the 
cubic splines. In the case of quasiuniform mesh we now get, as expected, ||?7— 7?/i|| = 0{h'^), \\u—Uh\\i = 0{h^), 
following the technique of the analogous proof in the previous section. Use of a uniform mesh and of relevant 
superconvergence results of Wahlbin, (3D] , allows proving a series of suitable superaccuracy estimates for the 
error of the interpolant and the error of the elliptic projection into the cubic spline subspaces. These results 
lead to error estimates such as ||?7 — %|| = 0{h^^^^ ^/In l/h), \\u — Uh\\ — 0{h'^^/lnl/h), that are consistent 
with the rates of convergence of the errors obtained from numerical experiments. 

In section 4 we turn to the analysis of fully discrete schemes. We consider only explicit time-stepping 
schemes in order to avoid the more costly implicit methods that require solving nonlinear systems of equations 
at every time step. Of course, with explicit methods there arises the issue of stability of the fully discrete 
schemes. We confine ourselves to a uniform mesh spatial discretization and consider three representative 
explicit Runge-Kutta schemes, namely, the Euler, the improved Euler, and the classical, four-stage, Runge- 
Kutta methods, whose orders of accuracy are 1, 2, and 4, respectively. We couple the Euler and improved 
Euler schemes with a piecewise linear spatial discretization and the fourth-order RK scheme with cubic 
splines. We show that the stability restrictions on the time step k needed by these schemes are of form 
k = 0{h^), k = 0{h^^^), and k < XqH for a constant Aq sufficiently small, for the schemes with Euler, 
improved Euler, and fourth-order RK temporal discretizations, respectively. Under these restrictions, we 
prove optimal-order in time error estimates for the fully discrete schemes; the spatial rates of convergence 
are these of the semidiscrete approximations. The evidence of numerical experiments is consistent with the 
theoretical results. 

In section 5 we analyze a nonstandard Galerkin semidiscrete approximation of (jCB|) and (jSCB|) that 
approximates 77 by piecewise linear continuous functions and u by piecewise quadratic functions defined 
on the same (uniform) mesh. (The method may be generalized using the analogous pairs of higher-order 
splines but here we confine ourselves to the low-order case.) A series of superaccuracy results for the error of 
the projection onto the space of piecewise linear continuous functions is developed, one of which makes 
use of a result of Demko, [17], on the exponential decay of the off-diagonal elements of the inverse of a 
certain tridiagonal matrix. With the aid of these results we show optimal-order error estimates e.g. of the 
type \\r] — r]h\\ = 0{h^), \\u — Uh\\ = 0{h^). Finally, in section 6 we make some remarks on the application 
of standard Galerkin methods on some simple first-order hyperbolic problems. Specifically, superaccuracy 
tools that were developed in the case of uniform mesh in previous sections are combined with further similar 
estimates proved in this section to provide the basis for extending the optimal-order error estimate of 
Dupont, [20], for the periodic initial value problem, to the case of two types of initial-boundary value problems 
for first-order hyperbolics. The situation is contrasted with what happens in the case of nonuniform meshes 
by means of numerical examples. 

In the paper we use the following notation: We let C'' — C'^[0, 1], k = 0, 1, 2, ... , denote the space of 
k times continuously differentiable functions on [0, 1] and define Cg = {0 G C''; (/)(0) = — 0}. We let 
for integer fc > H'', \\ ■ \\k denote the usual, L^-based Sobolev space of classes of functions on [0,1] and 
its associated norm. (In the case fc = 1 we use the equivalent norm defined by ||w||i — (||u|p -I- ■jH^'lP)"'^/^.) 
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The inner product and norm on = L^{0,1) we denote simply by || • ||, (•,•), respectively. The norms 
on = L°°(0, 1) and on = W^{0, 1) we denote by || • lU, || • |U,oo, respectively. We let be the 
polynomials of degree < r,and by (•,•), | • |, we denote the Euclidean inner product and norm on M.^ . 

2. Standard Galerkin semidiscretization with piecewise linear, continuous functions 



2.1. Semidiscretization on a quasiuniform mesh. Let = xi < X2 < ■ ■ ■ < xjv+i — 1 denote a 
quasiuniform partition of [0, 1] with h := maxi(xi+i — Xi), and let 5'? := {0 € : , i G Pi , 1 < J < 

N}, 5^ Q = {(/> e Sf^ , 0(0) — = 0}. Let Ihfi denote the interpolation operators, with respect to the 
partition {xj}, into the spaces 5^, Sf^ q respectively. Then, it is well known that there exists a constant C 
independent of h such that 

\\w- Ihw\\+h\\{w- IhwYW < Ch''\\w^'''>\\ , (2.1) 

for w € H'', k = 1,2, and that a similar estimate holds for Ih.ow if w S H'' n Hq. (In the sequel, C will 
denote a generic constant, independent of h.) Let a(-, •) denote the bilinear form 

a{^, X) := (V', X) + X') Vt/^, x e , 

and Rh : ^ S"^ q he the elliptic projection operator relative to a(-, •), defined by 



a 



It follows by standard estimates that for fc = 0, 1 

\\RhW~w\\k<Ch^~''\\w"\\ if wGH^DH^. (2.2) 

There also holds, [19], that 

\\RhW - w\\ao + h\\RhW - w||i,oo < C'/l^||w||2,oo , (2.3) 

provided that w G I^^.oo ^ 

As a consequence of the quasiuniformity of the mesh, the inverse inequalities 

llxlli <C//^"'llxll, (2.4) 

||xlloo<C/i-^/'||x||, (2.5) 

hold for any Sf^ (or x G Sf^ q), and so does the estimate, [19], 

\\Pv~v\\oo<Ch^\\v\\2,oo, for weM^2,oo^ (2.6) 

where P : — > S*^ is the L^-projection operator onto Sf^. 

The standard Galerkin semidiscretization on Sf^ of (jCBp is defined as follows: We seek 77^ : [0,r] — >■ 5*^, 
Uh : [0, T] ^ Sl^Q, such that for i e [0, T] 

(ryht, (/)) + {uhx,4>) + {(rihUh)x, </') = V0 S 5^ , 
a(w/it, X) + (?7fta:, x) + {uhUhx,x) = Vx e S'^^o - 

with initial conditions 

rih{0) ^ Pilo , Uh{0)^RhUo. 
Similarly, we define the analogous semidiscretization of (|SCBp . which is given for < t < T by 

[rjhtA) + {uhx,(l>) + \{{rihUh)x, 0) = V</> e Sl , 
a(u/jt,x) + {Vhx,x) + |(u/iM/ia;,x) + U'^h.Vhx,x) =0 Vx € 5^_o ! 
with 

Vh{0) = Pm, Uh{0) = RhUo. (2.10) 
Upon choice of a basis for S"^, it is seen that the semidiscrete problems (|2.7l) - (|2.8p and (|2.9p - (|2.10l) represent 
initial-value problems for systems of o.d.e's. Clearly, these systems have unique solutions at least locally in 
time. One conclusion of the next proposition is that they possess unique solutions up to t ~ T , where [0, T\ 
will denote henceforth the interval of existence of solutions of (jCBp or (jSCBp . 



(2.7) 



(2.9) 



Proposition 2.1. Let h be sufficiently small. Suppose that the solutions of (|Cj?[) . and (\SCB]\ . are such that 
r] e C(0,r; W^), u G C(0,T; Wl^nH^). Then, the semidiscrete problems ^1^, (HH and Hll), (flJ^ have 
unique solutions {r]h,Uh) for <t <T that satisfy 

max^||77(t)-%(t)|| <C/i, (2.11) 
ma,x \\u(t) -Uh(t)\\i < Ch. (2.12) 

0<t<T 

Proof. We first consider the case of the symmetric system (|SCB[) . Putting (/) = Vh and x = ''^/i in (|2.9p and 
adding the resuhing equations, we obtain the discrete analog of (ll.4[) . i.e. that 

h>m' + \\Mmi^\\vh{o)f + \\uhm\l (2.13) 

is valid in the temporal interval of existence of solutions of (I2.9p - (|2.10p . By standard o.d.e. theory we 
conclude that the system (|2.9p - (|2.10p possesses a unique solution in any finite temporal interval [0,T]. 

We now let p := 7] - Pt], Pi] - j]h, cr := u - RhU, £, := RhU - Uh. Using (|SCB)) and ([^ - (^111)) we 
obtain for < t < T, 

{etA) + {<^. + i.,4') + h{{m-iihUH)A)^Q y<l>esl, (2.i4) 

a{^t,X) + (Px + 6'rr,x) + |(ww:r - UhUhx , x) + HvV^ ~ VhVhx,x) = Vx e 5^,0 • (2.15) 

Note that 

f]u - rjhUh - u{p + e)-{p + e){a + 0+ ri{a + , 

VVx - VhVhx = [■np)x + {riO)x - {pO)x - PPx - OOx . 

By continuity, in view of (I2.10p . we conclude that there exists a maximal temporal instance th > Q such that 
\\0{t)\\oo < 1 for t <th. Suppose that th < T. Then, taking = 6* in (PH)) and using ^X^, (g^l), (1231), 
(|2.4[) . (|2.6[) . and integrating by parts, we have for < i < i/j 

hiPf - Hax,e) - {tx,o) -i[{{pu)x,e) + \{uxe,e) - ((pa),,^) 
- m)x,0) - Wxe^e) - \{s,x0,e) + ((w)x,0) + ((?70.,^] 
< Ik^llll^ll + + h\\Aoo\\Px\\\\e\\ + \\\ux\\oo\\p\\\\e\\ 

+ ill^i.llocll^^lP + ^llp.lllkllooll^^ll + \\\p\\\Wx\\oo\\0\\ + \\\p\\ooUx\m\ (2.16) 

+ cwpxmum + iik.iiooii^^ip + iiie.iiii^ii + ^hiiooiia.iineii 

+ ih.||oolk||||0|| + i||7?||oolie.lll|e|| + 
<c{h + u\\i + \\e\\)\\e\\, 

where C is independent of th. In addition, putting % = ^ in (j2.15l) wc similarly obtain, for < t < th 

+ \ [(w, ix) + {^Q. ix)-{pB,ix)- \{pix. p) - {eix.e)] 

< I|P||||C.|| + 11^1111^.11 + |hl|.oo|k||||e.|| + ilhi||oo||?||||e.|| 

+ ilkiiooiiaiie.li + ilkiiooik.iiiieii + ihiiooiiHiiie.ii 
+ ^iir?iiooiieiiiie.ii + iiipiiiie.ii + qiPiiiiiPiiiieiii + \Ux\m 



(2.17) 



<c[h' + m\i + \\e\\)m\ 



From and ((^Tf)) it is seen that for < i < th there holds 

im\ + m\i)<c{h + \\e\\ + m,), 

from which, by Gronwall's Lemma and ()2.10p . we conclude that 

\\e{t)\\ + \\i{t)\\i<ch, o<t<th, (2.18) 



where C is independent of th. Since by |16'||oo < Ch-^/^eW, if h is sufficiently small the maximality 

property of th is contradicted. Therefore we may take th = T, and (|2.1ip and (|2.12l) follow from (|2.18p . (|2.ip 
and 

Consider now (ICBI) and its semidiscrete analog (j2.7p - (|2.8l) . The invariance property (I2.13P no longer holds, 
and the i.v.p. (I2.7p - (|2.8p has a local unique solution. Using the same notation as in the case of (jSCBp . we 
consider the i.v.p. of finding 9{t) £ S^, ^{t) e Sf^ q ior t > 0, such that 

+ (a, + 0) + {[u{p + 9)-{p + 9){a + 0+ vi^ + 0]» = V<^ e 

ai^t,x) + iP.+0.,x)-{u{a + O-^i,x')-{^'J.+^Ux)=0 Vx e 5^,o > (2-19) 

e{o) = o, c(o) = o. 

Obviously, (|2.19p has a local unique solution. Let th E (0, T) be the maximal time instance for which this 
solution exists and satisfies ||6'(t)||oo < 1 for < t < i^- Then, as in the case of (ISCBp . we obtain again that 

\\e{t)\\ + umi < ch , o<t<th, 

where C is independent of th- We conclude that we may take th = T. If rjh — Prj — 6, Uh — Rhu — ^, it 
follows that {r]h,Uh) is a unique solution of (|2.7p - (|2.8p in [0,T] and that it satisfies the estimates (|2.1ip and 

2.2. Uniform mesh. We now turn to the case of uniform mesh. For integer iV > 2 we let /i = 1/N and 
Xi = {i — l)h, i = 1,2, . . . , N + 1. (We shall use the previously established notation for the finite-dimensional 
spaces S^, S'^g, and their associated interpolation and projection operators.) We let {(j^jjfj'i denote the 
basis of si satisfying (pjixi) = Sij,l<i,j<N+l. In the following lemma we collect results that will prove 
useful in the error estimates that follow. 

Lemma 2.1. (i) Let Gij — 1 < i,j < N +1. Then there exist positive constants ci and C2 such 

that 

ci/i|7|^ < (Gt, 7) < C2/i|7p V7 e R^+^ . 
(ii) Let b e 7 = G-^b, and = 7j</>j. Then 

\m<{c,h)-^/^\b\. 

(Hi) Let w £ . Then, there exists a constant Ci — Ci(||u''^^-'||oo) such that for any x £ [xi,Xij^i\ 

{w - Ihw){x) = ~\w"{x){x - Xi){xi+i - x)+ g{x) , Xi < x < Xi+i , 
where \\g\\oo + ^lls'lloo < Cih^ . 

Proof. The proofs of (i), (ii), (iii) are given in [50] in case the finite element subspace consists of continuous, 
piecewise linear, periodic functions on [0, 1]. It is straightforward to adapt them in the case of Sf^ at hand. □ 

The following superapproximation result, a consequence of cancellations due to the uniform mesh, is 
central for the sequel. 

Lemma 2.2. Let v E and w E . If e := v — IhV and t/j E is such that 

(^,(/.) = ((i«e)',0) Vcl^ESl, 
then WipW < Ch^l"^. If m addition w(0) = w{\) = 0, then \\tp\\ < CK^ . 

Proof. Let 6, := {{we)' = -(we,^^, 1 < i < iV + 1. Clearly, = ©(/i^) and |6Ar+i| = 0(/i2). Let 
2 < i < N . Then, from Lemma 2.1(iii) we have 

bi — " ^^^^ / w{x){x — Xi^i){xi — x)dx — ^ „^f''' / w{x){x — Xi){xiJ^.i — x)dx + 0{h'^) 



-v"(x ) r^^^ 

= '^2h I [w{x) — w{x — h)]{x — Xi){xi^i — x)dx + 0{h'^) 

J Xi 

u''(ex)(x - x,)(x,+i - x)dx + 0{h^) . 



Therefore, bi = 0{h^), 2 < i < N; consequently |6| = 0{h^) and = 0{h^/^) by Lemma 2.1(m). If 
w(0) = 0, Lemma 2.1{iii) gives 

= -{we, (t)[) = ^^^^ / w{x)x{h - + 0{h^) 



= / - x)dx + 0(/i3) = 0(/i3) 

Jo 

Similarly, if = 0, then bN+i = 0{h^). Hence, \b\ = 0{h^/'^), giving HV-H = 0(/i^) by Lemma 2.1(m). □ 

In the uniform mesh case at hand we consider again the semidiscretizations ()2.7p and p.9p of (jCB|) and 
(|SCB[) respectively, with the initial conditions for both systems now given by the interpolants of 770, uq: 

Vh{0) = IhVo, Uh{0) = h^oUQ . (2.20) 

The main result of this paragraph is 

Theorem 2.1. Let h = 1/N be sufficiently small. Suppose that the solutions of (jCi^P and f^SCB^ are such 
that ?7 G C{0,T;C^), rjt e C{0,T;C^), u, Ut e C(0,T;Cf^). Then, the semidiscrete problems 
and \2.9^ . \2.2U^ have unique solutions {rih,Uh) for <t <T that satisfy 

(i) max - 77/j(i)|| < Cft.^/^ max - it/j(t)||i < C/i , 

0<t<T 0<t<T 

(ii) max \\u(t) - Uh(t)\\ < Ch^ , max WuAt) - Uht(t)\\ < Ch^ . 

0<t<T 0<t<T 

Proof. We give the proof in detail in the case of (|SCB[) . where existence of solutions of the i.v.p. (I2.9p . (|2.20p 
for < t <T follows from (|2.13p . The proof in the case of (ICBI) folows from an argument analogous to that 
given in the proof of Proposition 2.1 and will be omitted. Let 

p-.^Tj- Ihi] , 9 -.^ IhTj - Tjh , cr := u - IhflU , ^ IhfiU - Uh • 

Note that 

rju - TjhUh = 7]a + u9 - 9^ + F , 



where 

In addition 

where 

and 

where 



F := rjS^ + up — pa — — 9a . 
VVx - VhVhx = -99^ + {t]9)^j; + Gx 



G -.^ T]p - p9 - \p^ , 



UUx - UhUhx = Fix , 



H := Mcr + < - cr^ - icT^ - if 



2"^ 2 

Then, from ((SCB|) and ([2:201) it follows for < t < T that 



{9t,(p) + (C., + (((1 + \Tl)a)xA) + \{{u9)x.4>) - \m)xA) 



(2.21) 
(2.22) 



a(6,x) + idx,x) + {Px,x) - \{09x,}() + \{{vd)x,x) + ^(G.,X) 
with 

0(0) = 0, C(0) = 0. (2.23) 
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(In the right-hand side of (|2.22l) we used the fact that for x G S',^ q a{(7t, x) — (^t, x): since (w'— (//i,ow)', x') = 
for V G Hq.) In order to show the estimates in (i), we put (j) = 9 and x = ^ in (|2.2ip and (12.221) . integrate 
by parts, and add the resuhing equations to get for < t < T 

Hm\' + lieil?) + (((1 + ^V)cr).,9) + + 

Using now the approximation properties of S'^ and S'^ q, integration by parts, and Lemmas 2.1 and 2.2 we 
see that 

(((1 + ir?)cT)., 0) I = I (P[((l + ir;)a).], 0) | < Ch'/'\\e\\ , 
\iiu9).,e)\^^\{u,e,0)\<C\\0\\\ 
\{F,,0)\ < \mh,0)\ + \{{upUO)\ + \i{pa),,e)\ + \{{pO.,0)\ + 

< c^dieiiiii^^ii + h^im + h^e\\ + hu\\i\\o\\ + h\\0\\^) 

<c{mi\m + h'\\e\\ + h\\e\\'), 

\iP.,o\<ch^ai, 
\{m.,o\<cmi\\oh 

\iG.,0\ <\im^.)\ + \ip0,^.)\ + 
<c{h^ch + h^e\m\i), 

<c{h'mi + uf + h'm\l + h'mi) 

\iPue)\ + \{auO\<cih'\\e\\ + h^^\\). 

Hence, from (|2.24p we conclude that for < t < T 

im\' + mi)<c{h'^^e\\ + h^ai) 

+c{\\er + m\i)<c[h^ + i\\er + \m)]- 

Hence, by Gronwall's Lemma we obtain 

mt)r + umi<c[h' + mo)r + \\m\\i]. 

from which, in view of (|2.20p . we get 

\\d\\ + U\\i<Ch^/'^ , 0<t<T, (2.25) 

and the estimates in (i) follow. In addition, from (j2.25p and the approximation properties of 5^ and 5*,^ q 
one may easily derive the following estimates of F, G, and H, that we note for further reference 

ii^ii<c^(iieii + /^'), 

||G||<C/^^ (2.26) 
\\H\\<Cim\+h'). 

We proceed now to prove the optimal-order error estimates in (ii). Equation (j2.22p may be written in 
the form 

6 = Rhv , (2.27) 

where v is the solution of the problem 

v-^v" ^-i9 + pU~^m,-^iG-^0'+3H),~at, x G [0, 1] , 

?;(0) = i;(l) = 0. ' 



Considering the weak form of (|2.28p in Hqj ^-nd using integration by parts in the right-hand side, (|2.25p and 
((2:26)) we see that 

\\v\\i<Ch^/^. (2.29) 
In order to derive a bound for \\v\\, let C G and V be the solution of the problem 

V-lV" = C, xe[0,i], 
V{Q) ^V{1) = 0. 

Then, by (P?^ 

{v, C) = a(«, V) = {9 + p, V) + i^, V) + i(G - + 3^^ V') - {at,V) . (2.31) 
From ([23T|) with = 1 we see that (Ot + Pt A) = 0, < t < T . Hence, 



(2.30) 



6 + p)dx = j p{x, 0)dx =. J = const . 

Therefore, if for {x,t) £ [0, 1] x [0,r] 

ry(x,t) := I {9{s,t)+ p{s,t))ds- xJ , (2.32) 



it follows that 7 e Hq and jx — & + P — J- Hence, (|2.3ip yields 

(v, C) = (7., V) + '^(m.X) + i(G - + 3^ _ + ^j^ V') - iat,V) 

= -(7, V") - i(7, r/T/" + ri^V) + i(G - + 'sH-r^p + r^J, V) - (a*, V) . 

Now, using (j2.25p . ()2.26p . the approximation and inverse properties of 5^, S^q, and elliptic regularity in 
([230]) we obtain 

1(1^,01 < ii7iir"ii + cm\\v"\\ + \\v'\\) + cih' + mww + ch^vw 
<cih^ + \h\\ + mmch 

and conclude that 

lkll<C(/.^ + ll7ll + lieil). (2.33) 

Let now W the solution of the problem 



(2.34) 



WiO) = W{1) = 0. 

Using (112Z1), (EUSl), elliptic regularity in (PTMl) . and the estimates (P?^ and gives 

(e,6) = a(W,^t) = aiW,Rhv) = a{RhW,Rhv) 

= a{RhV, RhW) = a{v, RhW) = a(w, RhW - T4^) + a{W, v) 
= a{v,RhW-W) + {^,v) < ChUWlHi + UWIM 

<c{h' + \h\\ + mmh 

from which there follows that 

i^lief <C(/i^ + ll7lP + lief), 0<t<T. (2.35) 

In order to obtain the required optimal-order estimate for ||^|| from (|2.35p we need a similar estimate for a 
suitable approximation of 7. For this purpose, observe that (|2.2ip yields, ioT < t <T and G 5^, 

(7xt, 0) + i((w^)., 0) = ~{wx,(t>) , (2.36) 

where w =^ ^ + {1 + ^rj)a - ^0^ + ^F; note that w G and that (j^:^ and (1^:^ give < C{U\\+h'^). 
Using integration by parts and the definition of 7 in (j2.36p yields for < i < T 

{^t,c^') + ^{u^.,(^')^~iw-^up+^uJ,<j)') "id^eSl. (2.37) 



Consider now the space Sf-^^ of discontinuous, piecewise constant functions relative to the partition {xj}. 
Given any ip e S^^ consider in (|2.37l) (j) £ S"^ such that (f>' = Hence we have for < t < T 

where if :=w-iup+iu J satisfies < C{\\^\\+h^). Taking now ip = Po7 in the above, where Pq is the 
L^-projection operator onto S^^, yields for < i < T 

hftWPolf + 5("7.' ^07) - -(^^,^07) • (2.38) 
Since ||Po7-7|| < Ch\\"f\\i, cf. e.g. (16.24) of [IS], and ||7||i < Ch^/^ by ^Tm . we have 
|(u7.,Po7)l = |(^i7.,i'o7- 7) - 5(^x7,7)1 < C{h\\j\\l + 
<C{h\h\\l + \\Poj~jr + \\Pojr) 

<cih^+\\Pojr). 

Hence, yields, for < i < T, that 

kTt\\Poir<cih'+\\Poir+\m- (2.39) 

Now, from (|2:35)) . since II7II < |lPo7 - 7II + II-P07II < Ch^^^ + \\Poj\\, we get for < t < T 

kiur<cih^+\\Po7\\'+ur)- (2.40) 

Adding (|2.39p and ()2.40p we finally obtain by Gronwall's Lemma and p.23p that 

\\Polf + U\\^<Ch\ 0<t<T. (2.41) 

Therefore, the first inequality of the conclusion {ii) of the theorem holds; in addition, by similar estimates 
with the ones used above, we also obtain 

\h\\<Ch^, 0<t<T. (2.42) 
Finally, we prove the second estimate of (ii). Let Z be the solution of the problem 

Z - \Z" = it, a; e [0, 1] , 
Z{Q) = Z{1) = 0. 

Then, by ((2:27)) 

Witf = a{Z,^t) = a{Z, Rhv) = a{RhZ, v) = a{RhZ -Z,v) + {^^v) . 
Hence, eUiptic regularity and ^(2^ . ((233| . ^(1A\\ . ((2^ give 

Ut\? <ch\\zu\v\\,^\M\M<ch^it\\, 

i-G. Il^tll < Ch^ , and the second estimate of {ii) follows. □ 

Remark 2.1. It is not hard to see that the conclusions of Theorem 2.1 hold if we take rjh{Q) as any approxi- 
mation of 770 in 5^ with optimal-order rate of convergence. However, Uh{0) has to be taken as Ih,oUo or 
RhUQ. 

Remark 2.2. The superaccuracy estimate = 0{h^^'^) of (|2.25p combined with (|2.3p and Sobolev's 

inequality yields the L°° estimate \\u — Uh\\oo — 0{h'^/^) for u. 

Remark 2.3. In the case of a quasiuniform mesh with h — maxi(xi+i — Xi), one may easily check that the 
analog of Lemma 2.1 still holds; however in the proof of Lemma 2.2 there is no cancellation anymore from 
adjacent intervals and the conclusion is just that |j'0|i < Ch. As a consequence, the techniques of the proof 
of Theorem 2.1 now yield that \\r] — r]h\\ — 0{h) and ||u — = 0{h), and, instead of the optimal-order 
estimates in [ii), that ||m — Uh\\ = 0{h^/'^) and \\ut — Uht\\ = 0{h^/'^). However, for the linearized problem 

1 „ (x,t)e [0,l]x [o,r], 

Ut+Tlx- gMxxt = , 

u{o,t) = 0, u{i,t) = 0, t e [o,r] , ^ ' ' 

r]{x,0) = r]o{x) , u{x,0) ^ uo{x) , x G [0, 1] , 
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the last two estimates may be improved to yield optimal order, i.e. to give \\u — = Oih?), 
0(h?). Our numerical experiments suggest that \\u — Uh\\ — Oih?) even in the nonlinear case, 
not been able to prove this thusfar. 



\\Ut - Uht W = 

but we have 



2.3. Numerical experiments. We first consider the case of approximations on a uniform mesh with h = 
1/N on [0,1]. Table \TJ\ shows the errors and the associated rates of convergence, in the L^-, L°°- and 

L-^-errors L°°-errors 



N 


V 


order 


u 


order 


40 


0.1894(-1) 




0.1749(-3) 




80 


0.6849(-2) 


1.467 


0.4259(-4) 


2.038 


120 


0.3761(-2) 


1.478 


0.1877(-4) 


2.021 


160 


0.2454(-2) 


1.484 


0.1051(-4) 


2.015 


200 


0.1761(-2) 


1.487 


0.6710(-5) 


2.011 


240 


0.1342(-2) 


1.490 


0.4652{-5) 


2.009 


280 


0.1066{-2) 


1.491 


0.3413{-5) 


2.008 


320 


0.8738(-3) 


1.492 


0.2611(-5) 


2.007 


360 


0.7328(-3) 


1.493 


0.2062(-5) 


2.006 


400 


0.6261(-3) 


1.494 


0.1669(-5) 


2.005 


440 


0.5430(-3) 


1.495 


0.1379(-5) 


2.005 


480 


0.4767(-3) 


1.495 


0.1158(-5) 


2.004 


520 


0.4230(-3) 


1.496 


0.9864(-6) 


2.004 



N 


n 


order 


u 


order 


40 


0.1126 




0.5916(-3) 




80 


0.5617(-1) 


1.004 


0.1509(-3) 


1.971 


120 


0.3741(-1) 


1.002 


0.6755(-4) 


1.983 


160 


0.2804(-l) 


1.002 


0.3813(-4) 


1.988 


200 


0.2243(-l) 


1.001 


0.2445(-4) 


1.991 


240 


0.1868(-1) 


1.001 


0.1701(-4) 


1.992 


280 


0.1601(-1) 


1.001 


0.1251{-4) 


1.994 


320 


0.1401(-1) 


1.001 


0.9582(-5) 


1.994 


360 


0.1245(-1) 


1.001 


0.7575(-5) 


1.995 


400 


0.1121(~1) 


1.001 


0.6139(-5) 


1.996 


440 


0.1019(-1) 


1.001 


0.5075(-5) 


1.996 


480 


0.9337(-2) 


1.001 


0.4266(-5) 


1.996 


520 


0.8618(-2) 


1.001 


0.3636(-5) 


1.997 



N 


V 


order 


u 


order 


40 


0.2440(+l) 




0.2449(-l) 




80 


0.1776(+1) 


0.459 


0.1192(-1) 


1.039 


120 


0.1466(+1) 


0.473 


0.7875(-2) 


1.022 


160 


0.1277(+1) 


0.480 


0.5880(-2) 


1.016 


200 


0.1146(+1) 


0.484 


0.4691(-2) 


1.012 


240 


0.1049(+1) 


0.487 


0.3902(-2) 


1.010 


280 


0.9725 


0.489 


0.3341(-2) 


1.008 


320 


0.9109 


0.490 


0.2920(-2) 


1.007 


360 


0.8596 


0.492 


0.2594(-2) 


1.006 


400 


0.8161 


0.493 


0.2333(-2) 


1.006 


440 


0.7787 


0.493 


0.2120(-2) 


1.005 


480 


0.7459 


0.494 


0.1942(-2) 


1.005 


520 


0.7170 


0.494 


0.1792(-2) 


1.004 



Table 2.1. Errors and orders of convergence. (jCB|) system, standard Galerkin semidis- 
cretization with piecewise linear, continuous functions on a uniform mesh. 



i/^-norms at T = 1, of the standard Galerkin approximation (using piecewise linear, continuous functions) to 
the (jCBp system with suitable right-hand so that its exact solution is given hy rj = exp(2f)(cos(7ra;) + a; + 2), 
u — exp(— a;i)a; sin(7ra;). The system was integrated up to T = 1 with the classical, four-stage, fourth- 
order explicit Runge-Kutta method (see Section 4) using a time step k = h/10. We checked that the 
temporal error of the discretization was very small compared to the spatial error, so that the errors and 
rates of convergence shown are essentially those of the semidiscrete problem (|2.7|) . (I2.20|) . The table suggests 
that the rates of convergence of rj and u approach 3/2 and 2, respectively, thus confirming the relevant 
estimates of Theorem 2.1. It also suggests that — = 0{h) (confirming the result of Theorem 2.1), and 
ll?? — ??/i||i = 0{h^^^), ||?7 — ?7;i||oo = 0{h), \\u~Uh\\oo — 0{h?). Table [2?2l shows the analogous errors and rates 
in the case of a nonhomogeneous (|SCBp system with the same exact solution 77 = exp(2i)(cos(7ra::)-|-a::-|-2), u = 
exjp[—xt)x sin(7ra;), integrated on a uniform spatial mesh up to T = 1 with the same time stepping procedure. 
The convergence rates are essentially the same. Table [273l a) shows some results from a computation of the 
exact solution 77 = exp(2t)(cos(7rx) -I- a; -I- 2), u — exp(a;i)(sin(7ra;) + — x"^) of a suitably nonhomogeneous 
version of the (jCB|) system with the standard Galerkin method with piecewise linear, continuous functions 
on the quasiuniform mesh on [0, 1] given by /i2i-i — 1.2Ax, ft.2i = 0.8Ax, 1 < * < -^/2, where hi — Xi+i — Xi 
and Ax — 1/N. Again the fully discrete problem was solved by the fourth-order accurate four-stage classical 
explicit Runge-kutta scheme with timestep k = Ax/10. We integrated the system up to T = 0.4 starting 
with the L^-projections of 770 and uq on the finite element subspaces. The temporal error was much smaller 
than the spatial error. Tabl d2.3y a) shows the L^-errors for and u and the associated rates of convergence 
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Z/^-errors 



I/°°-errors 



AT 


1 


order 


u 


order 


40 


0.7423(-2) 




0.3613(-3) 




80 


0.2678(-2) 


1.471 


0.8849(-4) 


2.030 


120 


0.1469(-2) 


1.481 


0.3907(-4) 


2.016 


160 


0.9579(-3) 


1.486 


0.2190(-4) 


2.011 


200 


0.6871(-3) 


1.489 


0.1399(-4) 


2.009 


240 


0.5236(-3) 


1.491 


0.9703(-5) 


2.007 


280 


0.4160{-3) 


1.492 


0.7122{-5) 


2.006 


320 


0.3408(-3) 


1.493 


0.5449(-5) 


2.005 


360 


0.2858(-3) 


1.494 


0.4303(-5) 


2.005 


400 


0.2441(-3) 


1.495 


0.3484(-5) 


2.004 


440 


0.2117(-3) 


1.495 


0.2878(-5) 


2.004 


480 


0.1859(-3) 


1.496 


0.2418(-5) 


2.003 


520 


0.1649(-3) 


1.496 


0.2060(-5) 


2.003 



N 




order 


u 


order 


40 


0.3491(-1) 




0.9655(-3) 




80 


0.1704(-1) 


1.034 


0.2537(-3) 


1.928 


120 


0.1125(-1) 


1.024 


0.1140(-3) 


1.972 


160 


0.8394(-2) 


1.018 


0.6442(-4) 


1.985 


200 


0.6694(-2) 


1.014 


0.4132(-4) 


1.991 


240 


0.5567(-2) 


1.012 


0.2873(-4) 


1.994 


280 


0.4764(-2) 


1.010 


0.2112{-4) 


1.995 


320 


0.4164(-2) 


1.008 


0.1618(-4) 


1.996 


360 


0.3698(-2) 


1.008 


0.1279(-4) 


1.997 


400 


0.3326(-2) 


1.007 


0.1036(-4) 


1.998 


440 


0.3021(-2) 


1.006 


0.8563(-5) 


1.998 


480 


0.2768(-2) 


1.006 


0.7197(-5) 


1.998 


520 


0.2554(-2) 


1.005 


0.6133(-5) 


1.999 



N 


»? 


order 


u 


order 


40 


0.1049(+1) 




0.2670(-l) 




80 


0.7383 


0.506 


0.1249(-1) 


1.096 


120 


0.6026 


0.501 


0.8131(-2) 


1.059 


160 


0.5219 


0.500 


0.6024(-2) 


1.042 


200 


0.4668 


0.500 


0.4784(-2) 


1.033 


240 


0.4262 


0.499 


0.3967(-2) 


1.027 


280 


0.3946 


0.499 


0.3388(-2) 


1.023 


320 


0.3692 


0.499 


0.2956(-2) 


1.020 


360 


0.3481 


0.500 


0.2622(-2) 


1.018 


400 


0.3302 


0.500 


0.2356(-2) 


1.016 


440 


0.3149 


0.500 


0.2139(-2) 


1.015 


480 


0.3015 


0.500 


0.1959(-2) 


1.013 



Table 2.2. Errors and orders of convergence. (ISCB|) system, standard Galerkin semidis- 
cretization with piecewise linear, continuous functions on a uniform mesh. 



N 


V 


order 


u 


order 


80 


0.1277(-1) 




0.7432(-4) 




160 


0.6383(-2) 


1.000 


0.1858(^4) 


2.000 


240 


0.4258(-2) 


0.999 


0.8259(-5) 


2.000 


320 


0.3194(-2) 


0.999 


0.4646(^5) 


2.000 


400 


0.2556(-2) 


0.999 


0.2973(-5) 


2.000 


480 


0.2131(-2) 


0.999 


0.2065(-5) 


2.000 


560 


0.1826(-2) 


0.999 


0.1517(-5) 


2.000 


640 


0.1598(-2) 


0.999 


0.1161(-5) 


2.000 


720 


0.1421(-2) 


0.999 


0.9177(-5) 


2.000 



(a) 



N 


V 


order 


u 


order 


40 


0.5852(-l) 




0.1693(-2) 




80 


0.2933(-l) 


0.997 


0.4271(-3) 


1.987 


120 


0.1942(-1) 


1.017 


0.1899(-3) 


2.000 


160 


0.1449(-1) 


1.019 


0.1068(-3) 


2.000 


200 


0.1155(-1) 


1.016 


0.6834(-4) 


2.001 


240 


0.9600(-2) 


1.014 


0.4745(-4) 


2.001 


280 


0.8213(-2) 


1.012 


0.3486(-4) 


2.001 


320 


0.7176(-2) 


1.011 


0.2669(-4) 


2.001 


360 


0.6371(-2) 


1.010 


0.2109(-4) 


2.000 


400 


0.5729(-2) 


1.009 


0.1708(-4) 


2.001 



(b) 



Table 2.3. L^-errors and orders of convergence. (jCB[) system, standard Galerkin semidis- 
cretization with piecewise linear, continuous functions on a quasiuniform mesh with '^^jl' = 
1.5 (a), = 150 (b) 



at T = 0.4. The data strongly suggest that \\r\ — r]h\\ = 0(\i) and ||m — — 0{h^), thus confirming the 
relevant theoretical result for rj (cf. Proposition 2.1 and Remark 2.2), and supporting the conjecture that 
the rate of convergence for u is actually equal to 2 even in the case of the nonlinear problem; recall 
that the technique of proof of Theorem 2.1 in the case of a quasiuniform mesh gives a pessimistic bound 
of 0{h^^^) for ||m — Uh\\, cf. Remark 2.2. These results are confirmed by the rates shown in Table I^TST b). 
which was obtained by integrating the same problem with the same method on the quasiuniform mesh 
on [0, 1] with /iiQi-9 — 0.02Aa;, /iioi-s = 0.05Aa;, /iioi-7 = O.OSAx, /iioi-e = 0.35Aa;, /iioi-5 — 0.5Ax, 
/iioi-4 = hioi-3 — Ax, /iioi-2 — hioi-i — 2Ax, hioi — 3Ax, 1 < i < N/10, and k — Aa;/10. 
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3. Standard Galerkin semidiscretization with cubic splines 

3.1. Semidiscretization on a quasiuniform mesh. Let = xi < 0:2 < . . . < a:jv+i = 1 denote a 
quasiuniform partition of [0, 1] with h :~ maxi{xi-i-i — Xi) and let 

St WeC': 0| e P3 , 1 < J < iV} , = Sf, : 0(0) = <^(1) = 0} , 

be the space of (the C^) cubic splines on [0, 1] relative to the partition {xj}, and the space of cubic splines 
that vanish at x = and at a: = 1. In this section we shall denote by //j : — > Sf^ the interpolation 
operator, with the properties that for any v e C^, {Ihv){xi) = v{xi), 1 < i < N + 1, {Ihv)'{xk) — v'{xk), 
k = 1, N + 1, and let Ih,Q ■ Cq — > Sf^ q be the analogous interpolant onto Sf^ q- It is well known that 

2 

J2h^\w- Ihw\\j <Ch''\\w^'''>\\ (3.1) 
3=0 

holds for any w E for fc = 2, 3, 4 and that a similar estimate holds for IhflW if w G n iJg. More 
generally, [28^, if 1 < fc < 4 and < j < fc we have that 

\\{w-x)^''^\\<Ch^-'-'\\w''^'^\\ if weH^ 



and 



mm 



min |l(«;-x)(^')||oo<C/i'=-^ ||wW|loo if w G C\ 



and that a similar estimate holds in Sf^ for w in H*' or in C'' that also vanishes at a; = and x — 1. Let 
a(-, •) denote again the bilinear form 

a(^, x) = (^, x) + x') , VV-, X e 5,t,o . 

and : iJ^ Sf^ q be the associated elliptic projection operator defined by 

a{RhV,x) ^ a{v,x) , Vx € 5^ q • 
It follows by standard estimates that if 1 < fc < 4 

\\RhW-w\\.i <Ch''-^\\w\\k j ^0,1, ii w € H'' nH^ , (3.2) 

and that, jl9) . 

\\RhW - wiloo + h\\RhW - u.||i^oo < C/i^||w||4,oo if w G W*^°^ n . (3.3) 

Similar estimates hold for the analogous elliptic projection into S^- In addition, as a consequence of the 
quasiuniformity of the mesh, the inverse inequalities 

\\xh<Ch-^P-"^\\x\U, 0<a<l3<2, 

\\x\\s,oo<Ch-^'+'/^'>\\xh 0<s<2, 

hold for any x & (or any x G S"^ g), and so does the estimate, [19], 

\\Pv - v\\^ < Ch^v\\4,oo , veW^'°^, (3.5) 

where P : ^ Sf^ is the L^-projection operator onto S^- As a consequence of the approximation and 
inverse properties of the cubic spline spaces we also have that P is stable in L°° and in H^, and that Rh is 
stable in and in nH^. 

We let the standard Galerkin semidiscretization on of (jCBp be defined as follows: We seek rjh ■ [0, T] — > 
St,Uh: [0, T] St 0, such that for t G [0, T] 



,4 (3-6) 



(ry/it, 0) + {uhx,(f>) + {{■rihUh)x, </<) = ^(f) e Sf^ 
a{uht, x) + {Vhx,x) + {uhUhx,x) = Vx G S'^^g 
with initial conditions 

??/.(0) = Pr?o, Uh{0) = RhUo. (3.7) 
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We also similarly define the analogous semidiscretization of (|SCB[) which is given for t E [0,T] by 
a{uht,x) + iWix,x) + i{uhUhx,X) + 2i'nh'nhxA) =0 Vx G , 

with 

Tlh{^) = Pm, UhiO) = RhUo. (3.9) 

The following result may be proved in a totally analogous manner to the analogous result of Proposition 2.1 
and, hence, its proof will be omitted. 

Proposition 3.1. Let h be sufficiently small. Suppose that the solutions of (j6'i?P and ^SCB\i are such that 
rj G C(0, T; W^), u G C(0, T; n Hi). Then, the semidiscrete problems ^^M, "■^'^ (E3), (E3 have 

unique solutions {rjh, Uh) for < t < T that satisfy 

max Mt)~r]hit)\\ < Ch^ , (3.10) 

0<t<T 

max \\u(t) - Uh(t)\\i < Ch^ . (3.11) 

0<t<T 

3.2. Uniform mesh. For integer N > 2 welet h — 1/N and Xi = {i — l)h, i = 1,2, . . . , N + 1. In this section 
we will denote by {(t)j}fJi^ the usual B-spline basis of S"^ defined by the restrictions on [0, 1] of the functions 
4>j{x) = $(| — (j — 2)), where $ is the cubic spline on R with respect to the partition {—2, —1, 0, 1, 2} with 
support [-2,2] and nodal values $(0) = 1, $(±1) = 1/4, $(±2) = 0. Thus, e.g. supp((?!)j) = [xj-3,Xj+i] 
and 4>j{xj-i) = 1 for 4: < j < N , etc. Our aim in this paragraph is to prove that the cubic spline standard 
Galerkin semidiscretizations of the two systems on this mesh satisfy the estimates ||?? — = 0{h^'^ y^\nl/h), 
\\u — Uh\\i = 0{h^), \\u — Uh\\ = 0{h'''y/lnl/h), \\ut — Uht\\ = 0{h'^^/lnTJh). For this purpose we shall first 
state and prove a series of auxiliary results. 

Our first lemma is a well known result, the cubic spline analog of Lemma 2.1. 

Lemma 3.1. (i) Let Gij — 1 < i,j < A'^ + 3. Then, there exist positive constants ci and ci such 

that 

ci/i|7p < (07,7) < C2/i|7p V7GR^+^ 
(ii) Let b G R^+3, 7 = G-^b, and = J2f=i Ij'^j- Then 

ll^ll <(ci/i)-i/^|6|. 

(Hi) Let w £ . Then, there exists a constant Ci = Ci(||w'^^^||oo) such that for any x G 

{w - Lhw){x) = j,w'^'^\x){x - Xi)'^{xt+i - x)^ + g{x) , x^ < x < x^+i , 
where \\g\\oo + h\\g'\\oo < Cih^ . 

Proof. The proofs of {i)-(iii) are given in |20) in the case of periodic cubic splines. It is straightforward to 
adapt them to the case of and Ih at hand. □ 

We next prove a superaccuracy estimate for P[(wey], where e is the error of the cubic spline interpolant 
of a sufficiently smooth function and w is a weight. This estimate is a consequence of cancellation effects 
due to the uniform mesh and may be viewed as the cubic spline analog of Lemma 2.2. 

Lemma 3.2. Let v G and w G . If e ~ v — IhV and 4> G is such that 

(V',0) = ((u;£)',0) V0G5^ 
then IIV'II < C/i^-S. If in addition w{0) = w{l) = 0, then < Ch^. 

Proof. Let bj — {{we)' , 0j) = —{we, (p'j), 1 < j < N + 3. In view of Lemma 3.1(m) it suffices to show that 
|6| = 0{h^). It is clear by (|3.ip and the properties of the basis functions (pj that bj — 0{h'^) for all j. We 



will prove that actually bj ^ OQv') for j = 4, 5, . . . , A^, thus establishing that \b\ ^ 0{h'^). Let A < j < N . 
Using Lemma i.liiii) and putting q{x) = x'^{h — a;)^/4! we have 

{we,(j)'j) —'^^ I w{x)e{x)(p'j{x)dx 

= w*'*^(xj_2) / w{x)q{x - Xj^3)(l)'j{x)dx + v''^\xj-i) / w{x)q{x - Xj^2)4''jix)dx 

+ v''^\xj^i) / w{x)q{x — Xj^i)if)'j{x)dx + v''^\xj) / w{x)q{x — Xj)(f)'j{x)dx + 0{h^) 
<y Xj-i J Xj 

V^''\xj-2)Jl+v'^^\x,.i){J2 + J3) + W^^'(X,)J4 + 0{h^) . 

Since u('')(xj_2) = t'^'^H^^j-i) ~ hv^^^tij) and u('')(a;j) = + hv^^'>{t2j) for some € (xj_2, 2:^-1) 

and G (^^j-ii^^j), we obtain 

{we,(t)'j) = v'-^\xj^i){Ji + J2 + J3 + Ja) + 0[h^) . (3.12) 
Suitable changes of variable in each of the four integrals yield 

Ji = w{x + Xj-3)q{x)(j)'^{x)dx , J2 — w{x + Xj-2)q{x)4'3{x)dx , 
Jo Jo 

nh ph 

J3 — w{x + Xj^i)q{x)(f>'2{x)dx , J4 = / w{x + Xj)q{x)(f)'i{x)dx . 
Jo Jo 

In addition, if Xj^ = {xj^2 + Xj-i)/2, we have for x e [0, h] 

W{X + Xj-3) = W{X + X.ji_i) - ^w'(rij^) , W{X + Xj_2) = w{x + Xj^) - ^w'{T2jx) , 
w{x + Xj-i) = w{x + Xj^) + |w'(T3ja;) , w{x + X-j) = w{x + Xj^) + ^w'{Tijx) , 

for some Tijx between x -\- Xj—4-\-i and x H- Xj^^ i = 1,2,3,4. Thus, 

Ji + J2 + J3 + Ji^ / w{x + Xjf,)q{x)[(t)\{x) + (t)'^{x) + 03(a;) + <l)'i{x)]dx + 0{h^) 
Jo 

h 



^w{x.j^) / q{x)[(j>'^(x) + (j)'2{x) + (j)'^{x) + (j)'^{x)]dx + 0{h') . 



The last integral is equal to zero, since 04(a;) = (f>i{h — x), 4>3{x) = (fii^h — x), and q(x) = q{h — x) for 
X £ [0,h]. Hence 

Ji + J2 + J3 + J4 = 0{h^) , 

and (|3.12p implies that bj = 0{h^). Thus, the first assertion of the lemma is verified. To prove the second 
assertion, suppose that w(0) ~ 0. Then bj — 0{h^) for j = 1, 2, 3. Indeed, using Lemma 3.1(m«), we have 

-5i = (u>£,0'i) = / w{x)q{x)(p[{x)dx + 0{h^) 

Jo 
l-h 

= v'^^\h) / xw' {tx)q{x)(t)'^{x)dx + 0(/i^) = 0{h^) . 
Jo 
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In addition, 

ph p2h 

—b2 — {we,(f>'2)— / w{x)e{x)(f>'2{x)dx + / w{x)e{x)(f>'2{x)dx 

Jo Jh 
ph p2h 

= v^'^\h) / w{x)q{x)(l)'2{x)dx + v'^'^'>{h) / w{x)q{h - x)(l)'2{x)dx + 0{h^) 

Jo Jh 
ph ph 

= v^^\h) / w{x)q{x)(f>'2{x)dx + v^'^'>{h) / w{x + h)q{x)(l>[{x)dx + Oih'^) 



Finally, 



J xw' {t^)q{x)(l)2{x)dx + v^'^\h) I {w{x + h) - w{x))q{x)(l)[{x)dx 
-6i + 0(/i'^) =0(/i5). 



f.2h p3h 

[we^cj)'^) — I ■w{x)e{x)(f>';^{x)dx + j w{x)e{x)(l)'^{x)dx + / w{x)e{x)(l)'^{x)dx 

J2h 



v^^\h) I w{x)q{x)(f>'^{x)dx + J w{x + h)q{x)(f>'2{x)dx 
w{x + 2h)q{x)<j)[{x)dx] + 0{h^) 



v'-^Hh) I xw' {t^)q{x)(l)'^{x)dx - bi - b2 +0{h^) = 0{h^) 



Similarly, if w{l) = we have bj — 0{h^) for j = + 1, + 2, iV + 3. Hence, if w vanishes at a; = and 
X = 1, |6| = 0{h^) and the second assertion of the lemma follows by Lemma 3.1(ii). □ 

We shall also derive a superaccuracy estimate for P[(u'e)'], where e is the error of the elliptic projection 
of a function v £ Cq and w is a weight. For this purpose, we first state two superconvergence results 
that follow from the analysis of Wahlbin in [30] , and which are valid for interior nodes whose distance from 
the endpoints of the interval is at least of 0{h\nl/h). 

Proposition 3.2. Suppose that w G Cg and let Vh = RhV he its elliptic projection onto 5*^0- Then the 
following hold: 

(i) There exists a constant C independent of h such that 

\{v ~ Vh)' {xi)\ < Ch'^\\v\\\Y^ provided dist{xi, dl) > Cihln , (3.13) 

where Ci is a sufficiently large constant independent of h. 
{a) If Xi, Xi^i are two adjacent nodes for which the second inequality in 1^3.13^ holds and e{x) = v{x) — 
Vh{x), we have 

e{x^+l) - e{xi) = 0(h'>)\\v\\wi^ . (3.14) 

Proof, (i) The estimate (|3.13l) follows from Corollary 1.6.2 of j50j (which is strictly valid when the elliptic 
projection is defined by ({i^, 0') — {v' , (j)') Vcp ^ q), combined with Theorem 1.3.1 of [30] that allows us to 
state the result for Vh = RhV defined as in section 2 of the paper at hand. 

(a) The cancellation property expressed by p. 141) is a consequence of the fact that e{x) = v{x) — Vh{x) 
may be represented in the form 

e{x) = Qi{x) + 0{h'')\\w\\w^{xi,xi+i) , x G [xi,Xi+i] , 

where Qi{x) is a polynomial of degree four such that Qi{xi) = Qi{xi-\-i). This representation follows from 
the remarks in Example 1.8.2 of [10] and by adapting the arguments of the proof in Section 1.8 of [50] 
(which are valid for Hermite cubics) to the case of the cubic splines at hand. Hence e{xi+i) — 
e[xi) = 0{h^)\\w\\w^^(xi,xi+i)i and p.l4p follows by the function values superconvergence estimate for elliptic 
projections given in Theorem 1.3.2 of [30]. □ 

In the next lemma we present some further formulas for the error v — R^v that will be used in the sequel. 
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Lemma 3.3. (i) Let v £ Cq, Vh — Rnv, e = v — Vh, and 1 < i < N. Then, there exists a constant 
C = C{\\w\\w^) such that for any x G [x^jX^+i], 

e{x) ^ ji{x) + jiv'''^'>{x){x - Xi)'^{xi+i ~ x)"^ + Si{x) , xe[xi,Xi+i], (3.15) 

where ||5i||oo < Ch^ and 7,; is the cubic Hermite polynomial interpolating the values of e and its derivative 
at the nodes Xi and Xi+i. 

(ii) In addition to the hypotheses in (i), suppose that Xi and Xi+i satisfy the second inequality in \3.1S^ . 
Then, there is a constant C = C{\\v\\\Y'^) such that for any x € [a;i,a;i+i] 

e{x) ^ e{xi) + jiv^'^\x){x - Xif{xi+i - x)"^ + 5i{x) , x€[xi,Xi+i], (3.16) 

where \\5i\\ 00 < Ch^ . 

Proof, (i) By the standard representation of the error of Hermite interpolation we have, since Vh G F3 in 
[xi,Xi+i], that for x G [xi,Xi+i] there holds 

e{x) - 7j = j,{x - Xi)'^{xi+i - x)'^v^'^'>{t^) 

for some t^ G {xi, Xi+i). Hence p. 151) follows from the mean- value theorem since v £ C^. 

(a) We let / = [0, 1], and if {a: e / : dist{x,dl) > Cih\nl/h}. In view of ((XT5|) it suffices to show that 

^i{x) = e{xi) + 0{h^)\\v\\y/^^ for Xi,Xi+i & K . 

Now 

7i(a;) = e{xi)Ai^i{x) + e' {xi)Bi^i{x) + e{xi+i)Ai^2{x) + e' {xi+i)Bi^2{x) , 

where 

Ao(x)= [1 + i?,i(x) ^ 
(x) = [1 + , B,2 (x) = . 

Hence, using p.l3p and p.l4p we see that 

ji{x) = e{xi)A^^i{x) + e{xi+i)A^,2{x) + 0(/i^)||w||iy5^ 
= eix^){A^,i{x) + Ai^2ix)) + 0{h^)\\v\\w^ 
^ e{x^) + 0{h^)\\v\\w5^ , 

which concludes the proof of the lemma. □ 

We are now ready to prove the required estimate for P[(we)'], where w is a weight. 
Lemma 3.4. Let w G Cq and w £ . If e ~ v — RhV and %jj e Sf^ is such that 

(^,0) = ((u;e)',0) V0e^^ 



then IIV'II < Ch^"^ ^/\aA/h. 

Proof. Let hj = {{we)' ,4ij), 1 < j < iV + 3. In view of Lemma 3.1(m) it suffices to show that \h\ = 
0{h^y/h)Ajh). It is clear by p.2p and the properties of the basis functions cjij that hj = 0{h'^). We will 
show that for 'most' of the the indices j it is true that hj = 0{h^). Let / = [0,1] and K — {x £ I : 
dist{x,dl) > Cih\nl/h}. Since Xi^{i- l)h, 1 < i < iV + 1, it is clear that ii Xi e K then 1 + Ci lnl//i < 
i < N + l-Cihil/h. Therefore X, e if if and only if z e M ^{u [1 + Ci In l//i, 1- Ci In l//i]}. 

Let ji ~ minAi and m = maxAi. We shall show that 

\bj\<Ch^, j = /i + 4,...,m-l. (3.17) 

Indeed, for such j we have 

—hj — {we,(f)'j) = / w{x)e{x)(j)'j{x)dx 

= 2^ je(xj^3+fc) / w{x)4)'j{x)dx+ I {e{x) - e{xj_3j^k))w{x)4)'j{x)dx\ . 

fc=0 -Ixj-i+k -Ixj^-i+k 
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Using now Lemma 3.3(iz) in each one of the second group of integrals of this expression and similar consid- 
erations as in the proof of Lemma 3.2, we see that 

-bj = e{xj-3)Ji + e{xj-2)J2 + e(xj„i) J3 + e{xj)Ji + 0{h^) , 

where Ji = J!^^^^^' w{x)(l)j{x)dx, 1 < i < 4. But since e{xi+i) — e{xi) + 0{h^) when Xi,Xi+i g K 
(by Proposition 3.2{ii)), and taking into account that the Ji are bounded independently of h and that 
Ji + J2 + J3 + J4 = 0{h) (by similar considerations as in the proof of Lemma 3.2), we finally obtain p.l7p . 
By the definitions of /i and m we now have 

M+3 m-1 Ar+3 

i=l J=M+4 j=m 

< C[(C2 + Ci In + /i^] < Ca/i^ In i , 

where C3 = CsdlwUvys , 7^*1)1 thus concluding the proof of the lemma. □ 

The main result of this section follows. 

Theorem 3.1. Let h — 1/N be sufficiently small. Suppose that the solutions of (jC-Bp and l\SCB^ are such 
that 77 e C(0,T;C5), rit € C(0,T;C4), u,Ut € C(0,T;C^). T/ien the semidiscrete problems ((IJ) and (fX^ 
wi</i initial conditions 

VhiO) = IhVo , Uh{0) = RhUQ , (3.18) 
/ia?;e unique solutions (rih,Uh) for <t <T that satisfy 

(i) max \\ri(t) - rih(t)\\ < Ch^-^\/]n^ , max \\u(t) - Uh(t)\\i < Ch^ , 

0<t<T V " 0<t<T 

(ii) max \\u(t) ~ Uh(t)\\ < Ch'^x/lni- , max WutU) — Uht(t)\\ < Ch'^ \/ln i- . 

0<t<T \ ti o<t<T V " 

Proof. Again, we give the proof in the case of (jSCBp . noting that the analogous proof for (jCBp follows as in 
Proposition 2.1. The technique of proof is basically the same as the one used in Theorem 2.1 and we shall 
give here the details that are different. We define now 

p := 77 - Iht] , 9 := Ihrj - rjh , a := u - RhU , ^ := RhU - uh , 

and, arguing as in the proof of Theorem 2.1, we have for Q <t <T 

him? + ml) + (((1 + h)<^).:0) + hiiue^S + h{F.,e) 

+ + + i(G.,0 + (3.19) 

where 

F := rj^ + up — pa — pS, — 6cr ^ 

G -.^ i]p - p9 - ip^ , 

H ■.= ua + u^-(7C~ icT^ _ 1^2 

Using the approximation properties of S"^ and 5^ q, integration by parts, and Lemmas 3.2 and 3.4 we have 

\{{{i + ^ijW),,e)\<ch^-''^/h^\\0\\, \{{ue),,e)\<c\\e\\\ 
\iF.,o)\ < ciimm + h^ew + \\er) , \ip,,o\ < ch^^w,, 

\{{voh,o\ < cmm , i(G.,e)i < cdieiiiii^ii + h'^uh) , 

\{H.,0\<c{Uf + h'U\\i), \{pue)\<ch'\\e\\. 

(Here and in the rest of the proof of this theorem we have omitted details that are similar to those in the 
analogous steps of the proof of Theorem 2.1.) Therefore, p. 191) gives for < t <T 



(ll^f + ||e||?)<C(/.^lni + ll^ll' + llell?): 



dt 
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from which, by Gronwall's Lemma and p.lSp we get 

\\0\\ + mh<Ch'-'^, (3.20) 

and the estimates (i) follow. 

In order to prove the estimates (m), we first note that that 

\\F\\<Cm\+h^), \\G\\<Ch\ <C(||e||+/^^). 
Noting that = Ruv, where v is the solution of the problem 

V - \v" = -{e + p)., - \{r^e)^ - \{G - + sh), , x e [o, i] , 

we see by (|3.20p that ||t;||i < Ch^^^^ ^ylnl/h. Arguing as in the proof of Theorem 2.1 we have now ||u|| < 
C(/i'^ + hll + ll^ll), where 7 = {0{s,t) + p{s,t))ds - x p{s,0)ds. Then there follows that 

^iur<cih'+\\Po7\\'+ur), 

where Pq is the L^-projection operator onto Sf^, the space of piecewise quadratic functions relative to the 
partition {xj}. We also have, in view of (I3.20p and the estimate ||Po7 ^ 7II < C^IItIIi, that 

'^i\\Poir<cihHn^^ + \\Po-,r+ur)- 

From the last two inequalities and Gronwall's Lemma we obtain 

\\Pol\\ + U\\<Ch^^/h^, 

from which the first estimate of {ii) follows. The second estimate also follows along the lines of the proof of 
Theorem 2.1. □ 

Remark 3.1. The results of the Theorem also hold if we take as ?7/i(0) any other approximation of ?7o in 
of optimal order of accuracy in L^. 

Remark 3.2. The superaccuracy estimate ||^|| 1 = 0{h^-^ ^/]nl/h) of (|3.20p . combined with p.3p and Sobolev's 
inequality yield the estimate \\u — u/i||oo — 0{h^'^ y/lnl/h). 

3.3. Numerical experiments. We considered the nonhomogeneous (ISCBp system and we discretized it 
on a uniform mesh with diminishing h = 1/N on [0,1] using cubic splines for the spatial discretization 
and fourth-order accurate explicit Runge Kutta time stepping, (cf. Section 4), with time step k = h/lO 
for which the temporal discretization error was negligible in comparison with the spatial error. We took 
a suitable right-hand side so that the exact solution of the system was 77 = exp(2i)(cos(7r2;) + x + 2), 
u = exp(a;t)(sin(7ra;)-|-x'^ — a;^). The errors and orders of convergence produced by this numerical experiments 
are shown in Table [XT] The rates are close to the theoretical predictions of Theorem 3.1. The table suggests 
that the rate of convergence for 77 is slightly less than 3.5, while that for u is essentially four. It 
further suggests that \\ri — rjfiWi = 0{h^'^), \\u — Uh\\i — 0{h^) (agreeing with the second estimate of (i) 
of Theorem 3.1), II77 — ?/;i||oo — 0{h^) and \\u — Uh||oo = 0{h'^). We also mention that the 1^^'°° orders of 
convergence (not shown here) were approximately equal to 2.2 for rj and 3 for u, and that the convergence 
rates from a similar experiment with (jCBp were practically the same. In Figure 3.1 we plot the quantity 
K := \\rj — ?7/i||/(^'^'^\/ln I//1) as a function oi N ~ l/h; here ||r/ — r/hW are the L^-errors from Table lOI We 
observe that k apparently approaches a constant close to 0.13 as N grows, which seems to be consistent with 
the presence of a slow- varying modulation of h^'^ as /i — ^ 0, such as (In We close this paragraph with 

a remark on the 'effect of the boundary' on the error estimates of Theorem 3.1. The proofs of Lemma 3.4 and 
Theorem 3.1 suggest that the accuracy of "tp in Lemma 3.4 and e.g. of \\r] — r]h\\ in Theorem 3.1 degenerates 
near the boundary of the interval. This is consistent with the results of the following numerical experiment. 
We integrated in time the ()SCBp system on [0, 1] with suitable right-hand side and initial conditions so that 
the wave elevation is given by the travelling Gaussian profile ri{x, t) = 0.5 exp[— 144(a; — 0.5 — 0.2t)^] and the 
velocity by u{x,t) = (){\/ri -f 1 — l)x(x — 1). (We use cubic splines in space and the explicit fourth-order RK 
scheme in time.) The support of the initial ?7-profile is effectively contained in the interval [0.3,0.5] and the 
wave moves to the right and starts crossing the boundary at a; = 1 at about t ~ 1.5 (see Figure 3.2). In 
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3.982 


440 


0.2323(-8) 


2.988 


0.1790(-10) 


3.985 


480 


0.1797(-8) 


2.950 


0.1256(-10) 


4.070 


520 


0.1412(-8) 


3.014 


0.9168(-11) 


3.932 



H^-erroTS 



N 


ri 


order 


u 


order 


40 


0.1298(-3) 




0.2014(^4) 




80 


0.2215(-4) 


2.550 


0.2540(-5) 


2.987 


120 


0.7923(-5) 


2.536 


0.7550(-6) 


2.992 


160 


0.3829(-5) 


2.528 


0.3190(-6) 


2.995 


200 


0.2182(-5) 


2.521 


0.1635(-6) 


2.996 


240 


0.1379(^5) 


2.516 


0.9467(-7) 


2.997 


280 


0.9363(-6) 


2.512 


0.5964(-7) 


2.997 


320 


0.6699(-6) 


2.508 


0.3997(-7) 


2.998 


360 


0.4987(-6) 


2.506 


0.2808(-7) 


2.998 


400 


0.3831(-6) 


2.503 


0.2047(-7) 


2.998 


440 


0.3018{-6) 


2.502 


0.1538(-7) 


2.998 


480 


0.2428(-6) 


2.500 


0.1185(-7) 


2.998 


520 


0.1988(-6) 


2.499 


0.9323(-8) 


2.999 



Table 3.1. Errors and orders of convergence. (jSCBp system, standard Galerkin semidis- 
cretization with cubic splines on a uniform mesh. 



0.18 
0.17 
0.16 
0.15 
0.14 
0.13 
0.12 
0.11 
0.1 
0.09 



Figure 3.1. k ■= \\r] - rih\\/h^ '^y^\nTJh as a function of = 1/h; the \\ri - r]h\\ are the 
errors from Table 

Table [5?^ we show the errors of 77, as = 1/h increases, at the temporal instances t = 0.5, 1.0, 1.5, 2.0 
and 2.5. The rates of convergence are practically equal to four up to t = 1.5 but as ij becomes nonzero at 
the boundary they fall to a value consistent with the first inequality of (i) of Theorem 3.1. 
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(a) »7h at t = 0.0 (b) Tjh at t = 0.5 (c) rjh at 4 = 1.5 

Figure 3.2. Travelling Gaussian T^-profile. Nonhoinogeneous (ISCBp system. 



time 


0.5 


1.0 


1.5 


2.0 


2.5 


N 


L^-error 


order 


L'' -error 


order 


-error 


order 


-error 


order 


-error 


order 


250 


1.3314(-08) 




1.0661(-08) 




1.3596(-08) 




1.5924(-08) 




1.9906(-08) 




500 


8.2780(-10) 


4.008 


6.6223(-10) 


4.009 


8.4585(-10) 


4.007 


1.0596(-09) 


3.910 


1.7594(-09) 


3.500 


750 


1.6334(-10) 


4.003 


1.3067(-10) 


4.003 


1.6706(-10) 


4.000 


2.2223(-10) 


3.852 


4.2637(-10) 


3.496 


1000 


5.1617(-11) 


4.004 


4.1350(-11) 


4.000 


5.2838(-ll) 


4.001 


7.4176(-11) 


3.814 


1.5595(-10) 


3.496 


1250 


2.1179{-11) 


3.992 


1.6922(-11) 


4.004 


2.1710(-11) 


3.986 


3.1966{-11) 


3.772 


7.1471(~11) 


3.497 


1500 


1.0287(-11) 


3.961 


8.1703(-12) 


3.994 


1.0554(-11) 


3.956 


1.6213(-11) 


3.724 


3.7803(-ll) 


3.493 



Table 3.2. L^-errors of rj and orders of convergence. Example of Figure [3?2] 



4. Fully discrete schemes 

In this section we turn to the study of some temporal discretizations of the o.d.e. systems represented 
by the standard Galerkin spatial discretizations of (jCBp or (jSCBp . such as (|2.2p or (|2.4|) . for example. We 
shall confine ourselves to explicit time stepping schemes in order to avoid the more costly implicit methods 
that require solving nonlinear systems of equations at each time step. Of course, with explicit methods there 
arises the issue of stability of the fully discrete schemes. We will not be exhaustive in our analysis but we 
will study as examples three simple, well known explicit Runge-Kutta temporal discretizations, that require, 
respectively, stability conditions of the type k — 0{h^), k = 0{h'^^^), and k < X^h for Aq sufficiently small, 
where k is the time step. 

4.1. The explicit Euler scheme. Let M be a positive integer, k — T/M denote the (uniform) time step, 
and put = nk, n = 0, 1, . . . ,M. We consider the standard Galerkin semidiscretizations with piecewise 
linear, continuous functions on a uniform spatial mesh on [0,1] with h = 1/N, given by the initial- value 
problems dSj]), ([2:201) and (I2J|), (|2:20l) in the case of the (|CB|) and the (jSCBp systems, respectively. We 
discretize the systems in time with the explicit Euler scheme. Hence, we seek for < n < M HJ^ £ S^, 
S g, approximations of the solution 77(2;, t"), w(a;, t") of the (jCB[) system, such that for < n < M — 1 



(4.1) 



{H]i+' - hj:, <t>) + kiuji,, 0) + kiiHi:u;:),),cf>) = o e ^^ 

a{Uj:+' - UJ:,x) + fc(i?L, X) + HUJ:UJ:,,x) = O Vx e si^ , 
with 

i/jj = hv'' , C/^" = Ih,ou" . (4.2) 
The analogous fully discrete approximation of the (jSCBp system is defined, for < rt < Af — 1, by 

iH;:+' - hj:, cf,) + kiuji,, 0) + ii{Hj:u;:),), 0) = o y^esi, 
a{uj:+' ~ uj:, x) + KHii^x) + fiKK,, x) + I{hj:hj:^, x) = o vx e si, , 



(4.3) 



with 



Let A:L^^ Si,, be defined for f e hy 



Hi hif , {/° = /„^ou° . (4.4) 
a{Af,x) = {f,x) Vx€52 , (4.5) 



i.e. as the discrete solution operator such that Wh ~ Af, where Wh is the standard Galerkin approximation 
in Si Q of the solution of the two-point bvp —-^w" + w — f,0<x< 1, w{0) — 'w{l) = 0. From (|4.5p we 
have immediately that 

P/lli <qi/ll-i, (4.6) 

where the || • norm is defined for f E by 

ll/ll-l = sup -j—r- . 

With this notation in place and letting as before P denote the projection operator onto S*^, we may 
rewrite the fully discrete scheme (|4.1I) . (|4.2p for (ICBI) as 

i/r 1 - hj: + kpuj:, + kPinj^uj:), = o , 

_ (/« + kAHJi, + kA{Uj:U]:j = 0, 
for < n < M - 1, with H° = hrf, = h.ou^ ■ Similarly, for (I5UB|) we have from gS]), gH) 

Hi+^ - ij« + fcPL/L + ip{hj:uj:), = , 

_ [/« + fcAi/^^ + fA{uj:uj:,) + lAiHjiHj:,) = o , 

for < n < A/ — 1, with H^jl = /ft77°, C^/^ = IhfiU^ ■ We will prove error estimates for the schemes (14.71) and 
(|4.8p by comparing iJ^' with Ih-qit^) and with Ihfiu{t^). For this purpose, it is useful to establish the 
following estimates of the truncation errors of the interpolants. (In the sequel we shall analyze mainly the 
approximation of the (ISCBp system; the analogous results for (jCBp follow as in Sections 2 and 3. Frequently, 
we shall suppress the x variable, denoting e.g. ri(-,t) by ri{t) etc.) 

Lemma 4.1. Suppose that the solution {r],u) of ^SCB^ is sufficiently smooth in [0,T]. Let H{t) = Ihrjit), 
U{t) = Ih,ou{t), and define ip{t) S Sf^, C(t) G S^g forO<t<T by 

Ht + PU., + ^PiHU)^^^, (4.9) 

Ut + AH^ + |A(f/[/,) + \A{HH^) = C • (4.10) 

Then 

WQWi <Ch\ \Mi<Ch\ 



hold for < t < T . An analogous result holds for (jCj?p . 

Proof. Subtracting the equations P{rit + Ux + ■^(uri)x) = and (|4.9I) . and putting p rj — Ihrj, a := u — Ih^u, 
we obtain 

P{pt + [(1 + h)'^]x + Kup). - i(p'T),) = ~V ■ 
Therefore, using the approximation properties of S'^ and S'^ g, and Lemma 2.2, we have 

\m < \\pt\\ + \\P[{i + ^v)'^U + U^M-W + IKHxIl < c{h^ + h""^ + h^ + h^) < ch^/^ . 

Similarly, since e.g. by Lemma 2.2 

||p[(i + ^^v)ctU\ < WPi^vtc^U + ||P[(i + ^,vWtU\ < ch^'^ , 

we have 

\\M < \\pt\\ + \\P[{i + hyUW + hWPi^pUW + \\{p<jU\\ < c{h^ + + h^ + h^) < ce/-^ . 

Note now that for any x G 5*^0, (|4.5I) and the fact that (crj, x') = yield 

a{A{ut - \utxx) - Ut,x) = {ut,x) + ^(utx,Xx) - a{Ut,x) = a{at,x) = {<^t,x) = a{Aat,x) ■ 

Hence 

A{ut - \utxx) - Ut ^ A(7t , 

which implies, in view of the second p.d.e. of (jSCBp that 

Adt + C/t + A{i]x + \uux + \mx) = . 



Subtracting now this equation from (|4.10p we see, after some algebra, that 

A{c7t +Px + f [(uCT)a; - crcr^] + ^[{vp)x - PP'A) = -C • (4-12) 
Therefore, using (|4.6p and the approximation properties of and Sf^ q, we obtain 

IICIIi < C^dktll-i + + Wiua-^aWU + Wivp-^p'U-i 

< C{\\at\\ + \\p\\ + \\ua\\ + + ||r;p|| + \\p^) < Ch" . 
Similarly, after differentiating (|4.12l) with respect to t, we see that 

Kt\\i<ch\ 

thus ending the proof. The same results hold for (jCBp of course. □ 

We now proceed to prove error estimates for the explicit Euler-Galerkin schemes (|4.7p and (|4.8I) . We 
begin by a consistency result. 

Lemma 4.2. Suppose that the solution {ri,u) of (|5'Cj?P is sufficiently smooth. Let := H{t"-) = Ihrjif''), 
[/" = [/(<") = Ihfiu{t"-), and define, for < n < M — 1, (5" and 62 by the equations 

Then 

max {\\S1\\ + WS'^Wi) < Ck{k + h^^^) . 

0<n<M-l 



The analogous result holds for (|Cj?P as well. 

Proof Using gH) and we have, with V" = V'l^"), C" = C(i"), that S[' = - i?" - /ci/^" + yfcV", 

(5J = - [/" - kUp + fcC". Hence, for < n < Af - 1, 

\\Si\\ + \\S^\\i < - - + /cllV-"!! + ||C/"+i - J7" - fc?7,"|li + k\\C\\i 

< C{k^ + kh^/^ + k^ + kh^) < Ck{k + h^l^) , 

by Taylor's theorem and (14. lip . □ 

Proposition 4.1. Suppose that the solutions {rj,u) of \SCB^ and (|Ci?P are sufficiently smooth on [0,T]. 
Then, if p, — k/h"^ , there is a constant C = C{p), which is an increasing continuous function of p, such that 

max II i/," - ,7(t")|| < C{k + e'^) , max \\U^ - < C{k + h) , (4.13) 

0<n<A/ 0<n<A'f 

where (HJ^, Uf^) satisfy ([73 '^^ case may be. 

Proof. Consider the case of (jSCBp . We use the notation of Lemmas 4.1 and 4.2 and put e" = i/" — Hf^, 
e" — ?7" — Uf^. Using the definition of 5", 82, and (|4.8p we obtain, after some straightforward computations, 
that for < n < M - 1 

= e" _ kP[el + i(i/"e"), - i(e"e"), + i(t/"e"),] + 5^ , (4.14) 
e"+i = e" - + f(?7"e"), - f (e"e^) + i(i/"e"), - ie"e:j] + <5^ . (4.15) 

It follows from that 



||^(C/"e"),||i < C||C/"e"|| < C||e"|| , (4.16) 

||A(iJ"e"),||i < C||i7"e"|| < C||£"|| . (4.17) 
Let now < n* < M — 1 be the maximal integer such that 

||£"||i + ||e"||i < 1, 0<n<n*. (4.18) 
Then, for < n < 71*, using (|4.6p . we have 

||P(e;' + i(7?"e"). - i(£"e").)|| < ||e2|| + i||(i/"e").|| + i|l(e"e").|| 

<C(||e"||i + ||e"||i||£"||i) <q|e"||i, ^ ' ' 

||A(e"e^)||i + \\A{e-el)\W < C(||(e")2|| + ||(£")2||) < C(||e"|| + ||e"||) . (4.20) 
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Using in (liTi)) and we have, for < n < 71* 

< ||^« _ |P(C/"e"),|| + CA;||e"|li + ll<5i II , (4.21) 
||e"+i||i < + Ck{\\e^\\ + ||e"|l) + M'ih ■ (4.22) 

Now 

lie" ^ |p(c/"e")^||2 ^ ||e"||2 _^ ^||p(c/«£n)^||2 _ ^({/"e")^) . (4.23) 
But, by inverse assumptions 

||P(t/"e"),|| < ||(C/"£").|| < 

In addition, 

|(£",P(C/"£"),)| = |(e",([/"£"),)| = i|(C/:e",£")| < C||e"f . 

Hence, (li:^ yields 

lie" - |P(C/"e")^|p < ||e"|p + C/i/fc||e"|P + Cfc||e"||2 , 

i.e. that 

||e"-|P(C/"e").||<(l + C^fc)||e'l, 
where is a generic polynomial in of degree one. Therefore (|4.21l) becomes 

||e"+'|| < ||e"|| + C^k{\\e^ + ||e"||i) + \\5'l\\ . (4.24) 

Hence, by (14.221) . (|4.24p . and Lemma 4.2 we obtain for < n < n* 

||e"+i|| + ||e"+i||i < ||e"|| + ||e"||i + C^fc(||e"|| + ||e"||i) + Ck{k + /i^/Z) . 

By Gronwall's Lemma we conclude that 

||e"|| + ||e"||i <C(Ai,T)(fc + /i3/2), 0<n<n* + l, (4.25) 

where C(/i, T) = exp((Co + Ci^)T). Taking h sufficiently small, we see from the maximality property ()4.18p 
of n* that we may take n* = M — 1, and the conclusion of the proposition follows from (|2.ip . The case of 
(|CB[) is entirely similar. □ 

Remark 4.1. The estimate (|4.25[) and Sobolev's inequality imply that ||e"||oo = 0{k + h?/"^). Therefore, 
max„ ||u(r) - U'^Woo - 0{k + /i3/2). 

Remark 4.2. Consider the linearized problem (|2.43p . In this case, the analogous fully discrete scheme is 

Hl+'-Hi: + kPUl, = Q, 

u]:+'-u;:+kAHj:,^o, 

for < n < M — 1, with = Inrf , = Ih,QU^ ■ Consequently, using the notation of the proof of 
Proposition 4.1, we now have for < n < Af — 1 the error equations 

£"+1 = e" - kPel + (5^ , 

from which there easily follows the estimate 

||e"|| + ||e"||i < C{k + h^'^), 0<n<M, 

and the conclusions of the analog of Proposition 4.1, without the stability restriction k — 0{h?). In other 
words, the linearized system is not stiff. This may also be verified by examining the spectrum of the spatial 
discretization operator of the semidiscrete linearized system: The latter may be written for < t < T in the 
form 

T]M + LhUh = , 
MhUht + Lhrjh = , 
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where the operators : Sf^ q 



Si, Lh : Sj-^ ^ Si Q, Mh '■ SIq^ SIq are defined by the equations 



e Si., 



e Si 



{Lh1p,4>) = {lpx,V) -viy JhfiT'P ^ 

(Lh<P, ^) = {(t>x, -0) V(/« e Sl, ip £ SI q 
{Mh^p,x) = ai^,x) yi^,X&Slo. 

Hence, the semidiscrete system may be written on S^ x S^q as 
where Wh = [r]h,Uh]'^, and 

Ah = 

Therefore, the spectrum of the spatial discretization operator coincides with that of the generalized eigenvalue 
problem 



AhWht 


+ BhWh = 


0, 


7 ■ 


, Bh = 


'0 Lh 


Mh 


Lh 



(4.26) 

Denoting the eigenfunctions as Vh = [Hh, UhY' , where Hh, Uh are elements of Sf^, Sf^ q, respectively, regarded 
as vector spaces over C, we have LhUh 



BhVh — —XhAhVh 

i]'^, where Hh, Uh a: 
-XhHh, LhHh = -XhMhUh, from which 



{LhUh, Hh) — —Xh\\Hh\\'^ , 
(LhHh,Uh) ^ -Xh{,MhUh,Uh) , 

where the L^ inner product for complex valued functions is defined as (/, 5) — f{x)g{x)dx. Now 



(4.27) 



Therefore, from (|4.27p 



{LhUh,Hh) — {Uhx,Hh) — —{Uh,Hhx) — —{LhHh,Uh) ■ 
-2tlm{Uhx,Hh) 



(4.28) 



\Hh\\^ + a{Uh,Uh) ■ 

We conclude that the spectrum of (|4.26p consists of purely imaginary eigenvalues (and is symmetric about 
the origin as —Xh is also an eigenvalue corresponding to the eigenvector Vh)- Moreover, it follows from (j4.28p 
that 

2||C/,,,||||i/,,|' 



|A/.|< 



\Hh\\^ + a{Uh,Uh) 



< C, 



i.e. that the spectrum is bounded by a constant C independent of h, implying that the linearized semidiscrete 
problem is not stiff. 

As we saw in Proposition 4.1, when the nonlinear (jCB[) or (jSCB[) system is discetized by the explicit Euler- 
Galerkin method, the mesh condition k — 0{h?) is sufficient for stability. In a numerical experiment, we 
solved the nonlinear, nonhomogeneous ()CB|) system using as exact solution 77(0;, t) = exp(2t)(cos(7ra;) + a; + 2), 
u{x,t) — exp(a;t)(sin(7ra;) + — x^) for x £ [0, 1], and discretizing by the explicit Euler-standard Galerkin 
method with piecewise linear functions and fixed N = 1/h = 400. When we integrated up to T = 1 using 
k = h? , we obtained an L^ error for 77 that was approximately equal to 2.2090(— 4). The accuracy degenerated 
when we took k = h°' with decreasing a < 2. For example, at the temporal instance closest to T = 1 the 
computations yielded the following magnitudes of the L^ errors of ry: 

k = h^-^ 3.8839(-4) 

k = h^-^ 1.1257(-3) 

k = h^-^ 3.6917(-3) 

k = h^-"^ overflow at about t = 0.8 . 
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4.2. The improved Euler method. We next study the temporal discretization of the initial-value prob- 
lems (|2.7p . (|2.20p and (|2.9p . (|2.20p by the explicit, second-order accurate 'improved Euler' scheme, that may 
be written in the case of the o.d.e. y' = f{t, y) in the two-step form 



j;«+l=y" + fc/(r + |,y"'l). 

Using notation analogous to that established in the previous paragraph, in the case of (|CBp we seek for 

< n < M e Si, Ull e SI q, and for < n < A/ - 1 H]^'^ e 5^, U]^'^ G 5^ q, such that 

HI'' + |PC/^, + Ip{HJ:UI:)^ = , 

U-'' + \AHl^ + \A{UlUi:) = , 

Ur^' K + kAHl-' + kA{U-''U--^) = , 
for < n < M — 1, with = IhTf , U^l — IhfiU^- In the case of (jSCBp the analogous equations are 

u--' -uj: + h^AHi, + f + ^,a{hj:hj:j = o , 

h;:+' rj: + kpu-;^ + lP{H-,'u-^% = , 

-UJ: + kAHl;^ + f A(C/r V^'i) lA{H",^'HrJ) = , 

for < n < A/ — 1, with = Ih'i]'^, Ujl = Ih.ovP ■ In order to study the consistency and convergence of the 
schemes, we let again i?" = = hvit"-), W = ;7(i") = //i,ou(t"), where (r?, u) is the solution of ([CB|) 

or (ISCBl) . and define, in the case of (|SUB|) . G S',^ x Sl ^ for < n < M - 1 by the equations 

In the case of (jCBp i/"'^, U"'' are defined analogously. Our consistency result is: 

Lemma 4.3. Suppose that the solution ir],u) of f^SCB^ is sufficiently smooth and let A — k/h. Define, for 
< n < M — 1, di , 62 by the equations 

gn ^ jjn+l _ ^ ;.^^n,l _^ . (4.33) 

Then, there exists a constant Ci = Ci(A), which is a polynomial of X of degree one, such that 

max (||(5"|| -I- ||<5?||i) < Cik{k^ + h^/'^) . 

0<n<A/-l 



The analogous result holds for (|Cj?P as well. 

Proof Let < n < M - 1. By (|OT|) . gH) and (|4l0l) we have 

From these expressions, after some algebra, we obtain 

where 

< ^Hl'Ul" - |(C/" 4- - + l^fDC" + xV'"C" • (4.35) 

Hence, by (|432|) . (|434l) . (j4J| . and the above we obtain 

gn ^ jjn+1 _ jjn _ j^jjn _ l^jju ^ j^^n ^ l^^n _ l^p^n ^ kp^n^ (4 3g) 
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Hn,ljjn,l ^ jjnjjn ^ k(^HH^)]^ + , (4.40) 



Now, (|4.35p . in view of (|4.1ip and the approximation and inverse properties of S"^ and Sf^ q, gives 

IWI|i<C(A:2|liJr||i||C/r||i + fc||^"||i(||{/"||i + fc||L/f||i) 
+ fc||C"||i(||if"||i + kWHl^W,) + fc2||V^"||i||C"l|i) 
< c(fc2 + kh-^h^/^{l + ck) + kh^{l + cfc) + k^h-^h^l^) 

Therefore, by Taylor's theorem and (j4.1ip we have 

\\Si\\ < c{k^ + kh^'^ + k^h^'^ + k^h^ + k{k^ + Xh^'^) < Cik{k^ + h^/^) , (4.37) 

where Ci is a constant that is a polynomial of A of degree one. (Such constants will be generically denoted 
by Ci in the sequel of this proof.) In order to estimate ||(52 ||i note that by (|4.34p 

UnAfjnA ^ ^n^n ^ |([/[/^)» + ^ (4.38) 

where 

:= ^U^Ur, - |((C/" + |C/r)C"), + ^CQ • 
By (|4.1ip and the approximation properties of Sf^ q we have 

\\w'^\\ <C{k^ + kh^). (4.39) 

Similarly, 
where 

By (14. lip and the approximation and inverse properties of Sf^ we have 

\Wi\\<C{k' + (4.41) 
By (11321), (|i35]) . (|i35)) . and dHO]), we see now that 

= ([/"+! _ [/" - fc[/r - ^Ul[) + fcC" + Cr - -T^^^' + f ^^2 + |^«^3 ■ 

Therefore, by Taylor's theorem, (lilD) . (gl]), (H^S)) . dUI]), 

11^2 111 < c(fc^ + fc/i^ + k^h^ + fc^/i^/^ + fc^ + k^h^ + k^ + Xkh^'^) < Cik{k^ + h^''^) , 
which, with (|4.37p . concludes the proof of the Lemma. The case of (jCBp is entirely analogous. □ 

For the stability and convergence of the fully discrete scheme it does not suffice to suppose that k — 0{h). 
The following result shows that the stronger condition k = 0(/i*/'^) is sufhcient. 

Proposition 4.2. Suppose that the solutions {rj^u) of \SCB^ and (jCi^P are sufficiently smooth on [0,T]. 
Then, if ^ — k/h"^/^, there is a constant C = C{fi), which is an increasing continuous function of fi, such 
that 

max \\HJ: - ,7(t")|| < dk' + h^^') , max ||[/," - uit-)\\, < dk' + h) . 

0<n<M a<n<M 

Proof We consider ((SCB|) . and put e" = i/" - HJ^, e" = [/" - U'^, 6*" = H'''^ - Hj^'\ and C = t^"'^ - ■ 
We will show that 

max (||e"|| + ||e"||i) < C{k^ + h^/'^) , (4.42) 

0<n<M 

from which the conclusion of the proposition follows. From (|4.30p . (|4.32p . (|4.33p we have, for < n < Af — 1 

^ _ j^p^n _ kp{^Hn,ljjn,l _ jjuA^n,!-^^ _^ gn ^ (4 43) 

and 

gn+i ^^n_ f^j^Qu _ _ Ul'^Ul^) - |A(i/"'ii/;''i - Hl-^Hl'^) + 5^ . (4.44) 

From Lemma 4.3 we have an estimate of ||(5"|| + ||(52l|i- Our goal is to obtain suitable estimates of the 
remaining terms of the right-hand sides of (j4.43p and (|4.44l) in terms of |je"|| + ||e"||i. To do this, note first 
that (|OT1) and ((^301) give 

r = e" - |Pe;' - |P(if"C/" - H'^U]i), . 
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61" = e" - |p« - |a;5* , (4.45) 

where 

p'l P(f/"£"), , (4.46) 

and 

:= -Pe^ + \P{H''e''). - iP(e"e"), . (4.47) 

Similarly, 

r - e" - - f A([/"t/: - - |yl(i/"i7: - • 

But C/"J7^' - C/,,C/;;'^ = (C/"e")^ - e"e2, - HIH]^^ = (i/"£")^ - e"eS. Therefore 

C = e" - ^uj^ , (4.48) 

where 

+ |A(J7"e"), - iA(e"e^) + iA(iJ"e"), - i^(e"e") . (4.49) 

In addition, by (I4.45P we have 

H^^ljn,! _ Hl'^U]l'^ = - |C/"p^ - f C/"lj'i' + , (4.50) 

where 

(f/">i - [/")6'" + ' - Q-^C ■ (4.51) 
From we conclude therefore that for < n < Af - 1 

e"+i = £" - + ^p'i - fc< + S'l , (4.52) 

where 

p^^ P(;7>?), , (4.53) 

and 

Lol Pe^ ~ yu^l, - |P(t/"c.r). + iPc^s". • (4-54) 

Finally, using the identities V^^^U^'^ -U^^'^U]^'^ = (L/"'iC")x-C"C and W^^^H^-^-Hl'^Hl;^ = (ff"'i6'")^- 
6*" 6*^', we obtain from (|4:44|) that for < n < M - 1 

e"+i = e" - kAeZ - ^A{{U'^'^C')x - CO - |^((-ff"'^6'")^ - r6»^) + (5J . (4.55) 

We now estimate the various terms in the right-hand sides of (|4.52p and (|4.55p . Let < n* < M — 1 be the 
maximal index for which 

||e"l!i + ||e"|!i < 1, 0<n<n*. (4.56) 
Then, by (|4.47p . the approximation properties of Sf^ and (|4.56p . we have, for < n < n* 

IKII < ||e"||i + C||i7"||i||e"||i + q|e"||i|je"||i < C||e"||i . (4.57) 
By (|4.49[) . the approximation properties of S*^, 5^ g, and (|4.6p there follows for < n < n* 

Ik2"lli < C(||e"|| + |le"|l + ||e"||i||e"|| + ||£"|| + ||£"||i||£"||) < C(||e"|| + ||e"||) . (4.58) 
Hence, by (|4.48l) . for < n < n* 

lirili<C(|l£"|l + |le"|li). (4.59) 
In addition, by (|4.45p . (|4.46p . (|4.57l) and the inverse assumptions we have for Q < n < n* 

Wl < ||e"|| +Cf ||£"|| +C||e"||i < (1 + CA)||£"|1 +qie"||i < CA(||e"|| + ||e"||i) , (4.60) 

where we have put again A = k/h; in the sequel, Cx will denote various constants that depend polynomially 
on A. Note also that in view of (|4.56p we have for < n < n*, from (|4.57p . (|4.46p and inverse inequalities 

\K\\i < ||e"||i + CA^IIp^'IIi + CfclKlli < ||e"||i + C\\\pn + CA||c^i'|| < Cx . (4.61) 

By (HSID, g311), (|ilT|) . glMl), (11351) and we have for < n < n* 

M\\ < c(||c/"'i - c/"||i||r|| + ||i?"'i||||rili + l|e"lllirili) 

< C{{k + h^)\\9"\\ + linii + lir lllirili) (4.62) 
<CA(||e"|| + ||e"||i). 
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Also, by dMni), (Hm, dmi), for < n < n* 

Mh < - c/"iiiiie"iii + iii^"''iiiiiriii + iie"iiiiiriii) 

<{k + h^)h-'Cxi\\e^ + ||e"||i) + + + CA(||e"|l + (4.63) 

<CA(||e"|| + ||e"||i). 

Hence, by (|4.54l) . (|4.58p . (|4.57p . (|4.62p and the inverse inequalities we have for < n < n* 

IKII < + Cfc(||£"|| + ||e"||) + CAIIe'li + C,i\\e-\\ + < C,{\\e^ + ||e"||i) . 

Therefore, in the right-hand side of (|4.52p we have for < rt < rt* in view of Lemma 4.3, 

II - fc< + ^^'11 < CAA:(||e"|| + ||e"||i) + C^Hk + h"'^) . (4.64) 

We embark now upon obtaining a sharp L^-estimate of the remaining term e" — ^p" + \-P2 in (|4.52p . We 
have 

Ik" - fpr + ^p^ir = + + Up2f - He\Pi) + f (£",P2) - ^{pi,P2) ■ 

Now, by 

(e",p?) = (£",([/"£"),) = i(C/,"£",e"). 

Also, by (0321), 

(£", P2) = -(£", U-pl) = -\\pn' + (C/^e", p?) , 

(pr,P2) = i(or,p'i)- 

We conclude that 

Ik" - + ^P^NP - lk"ll' + t^llP^f - |(f/,"£",£") + ^(c/re",Pi) - f^(c/>r,Pi) ■ (4.65) 
Now, by the approximation and inverse properties of q we have by (|4.46l) , (j4.53p 

|([/,"£",£")| < C|k"||2, 

I(c7>",p?)|<qk"||||pril <c/i-i|k"f , 

|(t/>?,Pr)l<C||pril'<C/i-2||e"f , 

||p^NI<c;i-2|k"||. 

Inserting these estimates in (j4.65p and recalling that p — k/h^/^ we are led to the inequality 
Ik" - |Pi + P2 IP < (1 + Ckp"^ + Ck + CkX + CA:A')|k"||2 < (1 + C^A:)||£"||2 , 

and, hence, 

lk"-|p'i+fp2ll<(l + C^fc)lk"lk (4.66) 
where, by we denote a constant depending polynomially on p; we have replaced Ca's by C^'s since 
A = h'^^^p < p. We finally obtain from gJH), (|4:B4|) and (j432l) . for < n < n*, that 

lk"+'|| < lk"|| + C^fc(lk"|| + ||e"||i) + C^k{k^ + h"/^) . (4.67) 

We now estimate the || • ||i norm of the various terms in the right-hand side of ()4.55p . For < n < n*, by 
gH), gJl, (!436l) . and (jiTSQl) we have 

\\A{{u--'n. - roiii < ciic/"^ir - i(r)'ii < cwch < c(ik"ii + ik-ik) . 

Similarly, for < n < n*, by (|4.60p . we have 

\\A{{H^'^0")^ - 0"e'^)h < c\\H"^^e" ~ i(r)2|| < CAiirii < CA(|k"|| + ||e"||i) . 

Therefore, using (|4.55p . (j4.60p and Lemma 4.3 we have for < n < n* 

||e"+i||i < ||e"||i +Cfc||r|| + CAfc(||£"|| + ||e"||i) +CAfc(fc2 + h"^/^) 
< ||e"||i + C^Afc(|k"|| + ||e"||i) + Cxk{k^ + h"/^) . 
Adding now (|4.67p and (|4.68p . we conclude for < n < n* that 

||£"+i|| + ||e"+i||i < (1 + C^fc)(||£"|| + ||e"||i) + C^fc(fc2 + /i3/2) . 
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(4.68) 



Using Gronwall's Lemma and taking h sufficiently small we conclude that n* may be taken equal to M — 1 
and there holds that 

+ ||e"||i < exp{Cf,T){k^ + h^^^) , 0<n<M, 
i.e. that (|4.42l) is valid; the conclusion of the proposition follows. □ 

Remark 4.3. Arguing as in Remark 4.1 we have again \\UJ^ — u(i")||oo = 0{k^ + h^/"^). 

Remark 4.4. If we repeat the numerical experiment in Remark 4.2 using the improved Euler method for 
time stepping, we obtain, for = \/h = 400, the evolution of the errors — ?7(t"')|| when k — h and 
k = h"^/^ shown in Table HTT] These results suggest that the condition k = h'^/^ is probably necessary as well 
for the stability of the scheme. If we repeat the experiment using k = h and the fourth-order explicit classical 
R-K scheme for time stepping (this scheme, as will be shown in the next paragraph, is stable when k = ft.,) 
we obtain errors approximately equal to those of the last column of Table HTT] For example at = 0.95 and 
t" = 1.0 the analogous errors are equal to 0.1810(— 3) and 0.1954(— 3), respectively. This implies that the 
errors of the last column of Table [4T] are essentially due to the spatial discretization. 



t" 


k = h 


t" 


fc = h^r-^ 


0.05 


0.9594(-5) 


0.05090 


0.9833(-5) 


0.1 


0.1857(-4) 


0.10179 


0.1890(-4) 


0.3 


0.5436(-4) 


0.30198 


0.5394(-4) 


0.5 


0.8538(-4) 


0.50217 


0.8393(-4) 


0.7 


0.1211(-3) 


0.70236 


0.1206(-3) 


0.8 


0.3949(-3) 


0.80075 


0.1448(-3) 


0.825 


0.2214{-2) 


0.82450 


0.1507(-3) 


0.85 


0.1445(-1) 


0.85165 


0.1572(-3) 


0.9 


0.8082 


0.90254 


0.1691(-3) 


0.95 


0.7706(+18) 


0.95005 


0.1811(-3) 


1.0 


over flow 


1.00026 


0.1963(-3) 



Table 4.1. L^-errors — the standard Galerkin method with piecewise linear 

continuous functions with h = 1/400 and improved Euler time stepping with k = h and 
k — h^/^] example of Remark 4.2. 



4.3. Fourth-order Runge-Kutta scheme with cubic splines. Our third example is time stepping with 
the classical, fourth-order accurate four-stage explicit Runge-Kutta scheme, written in the case of the o.d.e. 

y' = f{t, y) in the three-step form 

2/"''=2/" + |/(i" + iy"), 
2/"'' = 2/" + A;/(r + fc,2/"'2), 

We will couple this scheme with a space discretization that uses cubic splines on a uniform mesh on [0, 1]. 
We recall that the semidiscrete problem can be written, e.g. in the case of the (|SCB[) . in the form p.Sp or, 
equivalently, as follows. We seek (ry/i, Uh) G C^(0, T; Sf^ x Sf^ q) such that 

Vht + Puhx + hP{rihUh)x = , 

1 < t < T , 

Uht + Arjhx + ^A{uhUhx) + ^AiTjhrjhx) = , (4.69) 

r]h{0) = hvo , Uh{0)^RhUo, 

where P : ^ 5jJ is the L'^-projection operator onto S^, A : L'^ Sf^ q is defined as in (I4.5p . but now in 
S'^Q, Ih is the interpolant in Sf^, and Rh is introduced in Paragraph 3.1. The semidiscrete scheme for the 
(jCBp system is defined analogously. 

We start with an estimation of the truncation errors of the semidiscrete equations when applied to Ih'ni^) 
and Rhu{t) similar to that of Lemma 4.1. 
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Lemma 4.4. Suppose that the solution {rih,Uh) of \SCB^ is sufficiently smooth in [0,T]. Let Hit) = Ihrjit), 
U{t) = Rhu[t) and define ^ = ^p{t) e Sf,, C = C{t) e ^^.o. forO<t<T by 

Ht + PU., + \P{HU)^ = ^ , (4.70) 

Ut + AH,, + IA{UU,) + \A{HH,) - C • (4.71) 

Then, for j — 0, 1, 2, 3 

\\d{i^\\ < C/i^'Vlnl/Zi, \\d{C\\i<Ch'^, (4.72) 
/loZd /or < t < T . An analogous result holds for (|Cj?P as we/L 

Proof. Subtracting the equations 

Pijt + Pu, + \P{riu), = , 
Ht + PU, + \P{HU), = , 
we obtain, putting p :— rj — H , a := u — U 

Ppt + Pa, + \P{t^u - HU), = -ij . 

Since 

r]u - HU = r]{u - U) + U{r] - i?) = r]{u ~ U) ~ {u - U){r] - H) + u{t] ~ H) ^ f]a + up - pa , 
we have 

Ppt + PcTx + \P{ria), + \P{up), - \P{pa), = , 
and conclude from Lemmas 3.2 and 3.4 that 



\\di^\\<ch^-''^ha/h, J = 0,1, 2, 3, 

which is the first estimate in (|4.72p By the definition of A we have 

A{ut - \ut,x) = Ut ■ 

Subtracting now the equations 

A{ut - ^utxx) + At], + ^A{uu,) + ^Airjr],) = , 
Ut + AH, + UiUU,) + UiHH,) = C , 



we obtain 

But 

Hence 



Ap, + ^A{uu, - UU,) + \A{r^r^, - HH,) = -~C . 
uu, - UU, = {ua), - aa, , 7777^ - HH, = {rjp), - pp. 



Ap, + ^{A{ua), - A{aa,)) + ^{A{i]p), ~ A{pp,)) = -C, 
and fr'om (|4.6p and the approximation properties of S^, S'^q we conclude 

Hch<ch\ 

for j = 0, 1, 2, 3, thus proving the Lemma. The (jCBp case is entirely similar. □ 

We now define the fully discrete scheme. We consider only the case of (ISCBp . that of (ICBp being analogous. 
We let as usual Af be a positive integer, k — T/M and i" = nk, for n = 0, 1, . . . , M. For < n < M we 
seek {HJ^, UJl) e Sf^ x fi-^ „ approximations of 7/(t"), u(i"), and for < ti < M - 1 {Hl'\ Uj^'^) e Sf^ x S^^^, 
j = 1, 2, 3, such that for 6 < n < M - 1 

H--' HJ: + ka,P{U-'^~' + ^H-'^-'U-^-')^ =0, 

u-'^-uj:+ ka,A{H-r' + iK'^-'Kr' + iK'-'h]:^-') = o , 
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for j = 1,2,3 and 



(4.74) 



= 0. 



i=i 



where 



and 



-i?^, U;:^^^UJ:, ai-a2-l/2, 03-!, 6i=&4 = l/6, &2 = &3 = 1/3 , 



(4.75) 



i/°=77,.(0), f/^° = w^(0). 
We first investigate the consistency of this scheme. We define H, U as in Lemma 4.4 and put H" — H(t"), 
t/" = C/(t"). We let e S't W"-3 e S"^ q, for j = 0, 1, 2, 3 be defined for < n < M - 1 by the equations 

and 

Before proving our consistency estimates we introduce some notation that will simplify the computations. 
Let 

$ = C/+ii?C/, F ^ H^ + ^UU^ + ^HH^, $" = $(r), F'^ = F{e), 
and consider the functions 

- ^hx + I'-^h ^hx + 2-"/i ^/la; ' 

for i = 0, 1, 2, 3. With a slight abuse of notation we put for j = 0, 1, 2, 3 



We may then write (|4.70p and (I4.7ip as 
and (|i?75)) . dHHl), respectively, as 



Ut + AF^C, 



u;:'^-uj:+ka,AF--^-'^o, 



(4.76) 
(4.77) 

(4.78) 



for j — 1, 2, 3, and 



4 



(4.79) 



Finally, we may write the ralations (I4.75P in the form 



yn,j _fjn^ kajP<^i :>-^ = , 
^rn.J „ jjn ^ kttjAF''''-'-^ = , 



(4.80) 



for ] = 1, 2, 3. 
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Lemma 4.5. Suppose that the solution (jj^u) of \SCB^ is sufficiently smooth and let X = k/h. Define, for 
< n < M — I, 6i and 62 by the equations 

4 

gn ^ jjn+1 _Jjn^ E ^J*"'^"^] ' (4.81) 
J = l 
4 

^ jjn+1 _ jjn ^ kA]^ h^F''^^-^] . (4.82) 

Then, there exists a constant C\ which is a polynomial of X with nonnegative coefficients, such that 



max \\S['\\ <Cxk{k^ + h^-^^\nl/h) , max \\5^\\i < Cxk{k^ + h'^) . 

0<n<Af— 1 0<n<Af— 1 

Proof. We need to find formulas for the functions T^"'-' , W'"''^ , for j = 1, 2, 3 up to some remainder terms. 
From the first equation of (|4.77p we have 

whence from (|4.73p 

=F" + aifciyr-aifcV'"- (4.83) 
Also, from the second equation of (|4.80p we obtain 

and so from (f4l7\\ 

W""'^ ^ + a^kUl" - aikC . (4.84) 

Thus 

ynAy^rn,! ^ ^njjn _^ aik{HU)'i + ojk^ HJ^UJ^ + < , (4.85) 

where, in view of (I4.72p and inverse properties 



IKIli <CA/i3-Vlnl//i. 

In addition, from (|4.84p . (|4.85l) we obtain 

y^rn.l lynSy^n,! ^ Jjn ^ ^H^'U" + aik{Ul' + \{HU)'^) + ^k^Hl'Ul" + , 

i.e. 

$"4 = $" + aik<i>'l + ^k^H^Uff + , (4.86) 

where 

Ikzlll <CA/i'-'Vl^- 

Since now 

we have from (|4.86p that 

yn,2 ^Hn_ „2fcP$^ - aia2fc'P$^t - ^k^PiH^Un, - a2kP{v^), , 
from which, in view of (|4.73p . (|4.72p . and inverse properties there follows that 

with 

Ui\\<Cxh\ 

and 

In addition, from (|4.84p 

= U"U; + aik{UU^)'i + ajk^Ul'Ull + w'l , 

where, in view of (|4.72p 

Ikill < Ckh\ 
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Similarly, from (|4.80p there follows 

where, in view of (|4.72p and inverse properties 

ik^ii-i <CA/l^ 

and 

Therefore, from the above relations and (|4.83p we have 
i.e. 

pn,l ^ pn _^ ^^^pn _^ ^2^2(3f/«[/n _^ ^HJ^HH) + , (4.88) 

where 

< = -ki^^, + f < + iu;^ . 

Since 

there follows from (|4.88l) . (I4.77p . and previous estimates that 

W"^^ = U" + a^kUl' + a,a2eui\ - aja^k' [|A(C/,"C/rj + 1^(17," ffj^)] + Ci" , (4.89) 

where 

||Cl"l|l<CAfc/i^ 

We now have from (|4.87p . (|4.89p . the fact that ai — a2, and similar estimates to those used above, that 
^«,2^„.2 ^ jjnjjn ^ ^^^(iJC/)^ + aia^k^H'^U^ + Hl'Ul' + i?t"C/") + amlk'^iH^'Uli + H^^tU^) 

where 

ll^alli < Cxh'''^^^ 

follows from similar, as above, estimates and from the stability of P in H^. Rewriting the previous equation 
as 

p.n,2^n,2 ^ jjnjjn _^ ^^^(ifC/)^" + aia2k^{HU)1^ - aia2k^ H'^Ul' + aialk^{HtUt)'i 

we have, by the definition of and (|4.89p 

^na = $" + a2fc$r + aia2fc^«'rt - ^k^H^U^ + ^k^{HtUt)1 

~ ala^e [iA(C/rC//;) + \A{H^H:,) + \H-A{U:U:,) + \H-A{H^H^^) (4.90) 

+ i{/"p(i/rt^r).]+cr + K- 

Since now 

we have from and (|i?7S)) that 

yn,3 ^ ^« _^ _^ a2A:2^;j + aia2k^H^tt + ^k^PiHl'Ul")^ + V'2 , (4-91) 

where 

\m\<cxh\ 

and (using the fact that ||A/||2 < C||/|| for / e i^), 

llV^alli <^^A/i'-'v/h^. 
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In addition, from (|4.89p 
and so 

where < Ca/i^. In addition, from (j4:87| 



2 

iKii-i <CA/l^ 



Hence, 
where 

and 

From the definition of F"'^^ (gST]), (jilOS]) and (jilOS]) we have 



where 

Hence, since 
we obtain 

where 



+ aia2fc^ [lAiUl^Ull) + '^A{H-Hll)] + Q , 



'5 ' 



IIC2II1 < Cxkh'^ + Ck^. 

From (|4JT|) . (jilOS]) we get 

+ aia2A;3[|i7"A(;7f C/^J + \H-A{H^H^^) + \U-P[H^U^),,] + . 
where 

Hence, from the definition of (j4.95l) and the last relation, we see that 

$«-3 = $" + + a2k^^u + aioafc^^rtt " ^k^{HtUt)1 

+ aia2fc3[|^([/rf/rj + i^(^r^rj + \h-a{u-u-^ 

+ lH-A{H^Hr,) + lU-P{H^urh] + Q + i«5" • 
We are now ready to estimate 5" . By (|4.90p and (|4.96p we have 

63$"^' + 64$""' = {hz + 64)*" + (^302 + 64)fc$r + (aia2&3 + a264)fc'*?t 

+ a^a2b^e^i^ - ^en^u^ + 63(cr + ^^^3") + HQ2 + h: 
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In addition, in view of (|4.86p 

+ 62$"^^ = (&i + 62)*" + 0162^:$^ + ^PHl'Ul' + b2V^ . 

Hence, 

4 

where 

V- = b^v- + 63(Cr + I'^a ) + h{Q + H) ■ 
Therefore, by (|4.8ip and (|4.76p we obtain 

where, by (14721) 



Ik"!! < CAfc/i^ Vlnl//i. 

Since, in addition 

IKIli <CA/i'-'v/ha7)^, 

we conclude that 



(4.97) 



WKW < c\k{k^ + h''-^^/h[T/h). 

To estimate 62 we must also compute F"'"^. To this end, using (j4.95p we have after some algebra 

+ aia2fc'[|C/M([/r[/;^) + iC/"A(i?f + < , 

where 

M\\<Cxh\ 

Similarly, from (|4.9ip . we have 

yn,3yn^3 ^ ffu jjn ^ k(^HH^)^ + a2k\HH,)^, + fli a2 fc' (^ ff. 

- aia2k'{H^Hl^)^ + ^k^[H" P{H'^Un,], + , 

where 

||f7 ll-i < Ca/i" , < Cxh^-''.yh^. 

Hence, from the definition of F"'^, (|4.9ip . (|4.97p and the last relation, we obtain 

+ a,a2k%P{H-un. + lU-A{uru-) + |[/M(i?rfft" ) (4.98) 

+ ii?"F(ifrc^,")j^+<, 

where 

We may now estimate J". By (|4.94p and (|4.98p we obtain 

63i^"'' + 64^^"'' = (&3 + bi)F" + (02^3 + bi)kF^ + (0102^3 + a2&4)fc'i^t" + aia264fc'i^t"t 
- a^a2b,k\IUl^Ull, + \H^Hll) + 63^ + ^4^^? • 

Also, in view of (|4.88l) 
Therefore 



where 

= 62W3 + 63U'5 + 64W? 
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Hence, from (|4.82p and (|4.77p we obtain 

5- = - [/« - fcC/« - ^U- - ^U-, - ^U-,, + kAw- + < , 

where, by (|?:7^ 

\\w^\\i < Ckh'^. 

Since we also have, by previous estimates, 

\Wi\\-i<Cxh\ 

we finally conclude that 

\\5^\\i<Cxk{k^ + h^), 

and the proof of the Lemma 4.5 is completed. □ 

We show stability and convergence of the fully discrete scheme in the course of the proof of the following 
result. 

Proposition 4.3. Suppose that the solution (rj^u) of ^SCB]\ is sufficiently smooth on [0,r]. Let A = k/h 
and (H'P^^U'fli) be the solution of l \4.73^ - (4-'^4h Then, there exists a positive constant Aq and a constant C 
independent of k and h, such that for A < Aq, 



max - Hi:\\ < C{k' + h'-'y/hO/h) , max \\u{t-) - U^li < C{k^ + h^) . 

S ^ -^'J < n < Af 

Proof. It suffices to show that 

max - H]l\\ + llf/" - UrAi) < C(k^ + h^-^^^^) . 

0<n<M 

Let 

e" = - HJi , e" = [/" - UJi , e"^^ = F"'^' - Hj^'^ , e"'^' = - U]^'^ , j:^ 1,2,3. 
Then, by (|475l) and (|4J8)) . we have for j = 1, 2, 3 

= e" - /cajP($"^J-i - , (4-99) 

e"^J = e" - fcajyl(F"^J-i - K''^^^) , (4.100) 

and by g^g]), and 

4 

^n+i ^ ^« _ fc^6^p($"'J-i _ + > (4-101) 

4 

gn+1 ^ g" _ fc^ 6j A(P"'^-i - + . (4.102) 

i=i 

We embark upon the proof by deriving suitable estimates for e"'^ and e"^^. To this end, putting $}J := 
we have 



and since 

we will have 
i.e. 

where 
Therefore, 



-^l^ e" + i(i/"e" - £"e") + iC/"e" , 
7^ = e»+ i(77«e"-e"e"). 

P(ci>" - $^), = + ^Pi , (4.103) 
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where 

Hence, from (I4.99P 



7r = Pii 



(4.104) 
(4.105) 



Let < < M — 1 be the maximal integer for which 

lk"||i + ||e"||i < 1, 0<n<n*. 

Then, for n < n*, 

||7rili<C(||e"|| + ||e"||i), 

and, in view of (|4.104p . 

hril <C(lk"ll + l|e"||i). 

Therefore, by (giMl), for < n < n*, 

lk"'i||<CA(||e"|| + ||e"|U), 

and also 

||e"^^||i<l + CA||7ril+CA||pril<CA, 
where, as before, C\ denotes a constant that is a polynomial of A with positive coefficients. 



(4.106) 
(4.107) 
(4.108) 



In order to estimate e"' we note, setting = ' , that 

h 



Since 



= (C/"e"). - e"eS , 



and 

we have 

and therefore, for n < n 



Tjn Tjn Tjri Tjn (Tjn^n \ 

t -t,^ ^e.^ + ^(U e )x-2^ e^ + ^{H e )x - 2^ ^x 



\\F--Fj:\U<C{\\e^ + \\e^). (4.109) 
Recalling the expression for e"'^ in (j4.10ip we have for n < n* 

||e"-i||i<C(||e"|| + ||e"|U). (4.110) 

We now turn to deriving estimates for e"'^ and e"^'^. With this aim in mind, we note for later reference that 
it follows from the definition of 1/"^', M^"'^ i^"^' that for < n < Af - 1, j = 0, 1, 2, 3 

||^"'^||i + ||T4^"'^||i<Ca, < Ca . 

Since now 

_ Hl^^U^^^ = y"'i(P^"'i - [/;;'i) - (W^"'i - C/,"'')(F"'i - i?^'') + VK"'i(l^"'i - il,"'') 
using (|4.80p we have by (|4.105p 



and 



and so 
where 



$".1 _ - ^C/"p? +72 : 



7^ = e"^i + iy"'ie"'i - ie"'ie"'i - ^?7"7]" - ^(Ai^")e"'i , 
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(4.111) 



Recalling the definition of p" in (14.1041) we obtain 
where 

By (|4.11ip and (|4.108p we deduce that 

||72lli<CA||e"'i||i + Cfc||7ri|i + CAA||e"'i||i. 
Hence, by (|4.110p and (|4.106l) . (|4.107p we have 

||72"lli<C^A(|k"|l + |le"||i), 

and therefore, by (j4.113p we see that 

||72"ll<CA(|k"|l + |le"||i), 
holds for n<n*. Now, by (|TSO)) and (HTSl) 

gn.2 ^ ^„,2 _ jjn.2 ^ _ ^^^p^^n,! _ ^ 

and from (I4.112P 
Thus, 

Ik"''||<||e"||+Cfc||p5'||+Cfc2||p^||+Cfc|l72"||, 
and using (|4.113p . (|4.104p and (|4.114p we conclude that for n < n* 

\\s^'^'\\<C,i\\e"\\ + m), 

and 

||£"^'||i<Ca. 

We next consider e"'^. Since 
and 

we obtain by (|4.108l) , for n < n* 

\\F--'~F-''\\^,<Cx{\\e--'\\ + \\e--%), 
giving, in view of (|4.107p and (|4.110p for n < n* 

Therefore, by (|4.100p we conclude that for n < n* 



„n.2l 



e 



i<C,{\\e-\\ + \\e-\U) 



To obtain analogous bounds for e"'^, e"'^ note that 
and 

^«.2p^«.2 _ ^n,2^n,2 ^ -^^„,2^^n,2 _ jjn,2sj _ ^p^„,2 _ jjn,2.|^^„,2 _ ^n,2sj ^ p^„,2(-^n,2 
— yn,2gn,2 _ ^n,2^n,2 _|_ '^n.2^n,2 

By (ITSO)) and (|4.115p we see that 
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Therefore 
where 

So, by KTU^ and (|I1T^ we have 



2 2 2 ^ ' ' 



where 

Now, by (|4.120p and (|4.117p we have for 71 < n* 

||7?lli<CA||e"^2||i + Cfc||7^N|i + CAfc||£"-2|li. 
Hence, by (|4TT9)) . (|4ll4| and (|4Tl6)) . 

||73"lli<CA(||e"|| + ||e"||i), 

and thus 

||73ll<CA(||e"|i + lle"||i), 

for n < 71*. Now 



(4.120) 

(4.121) 
(4.122) 



„n,3 _n k n i 02 h n aia2k n 



Pi ^ ^P2 



-P? - ^7? ■ 



and from ()4.12ip 
Therefore, 

and by (|4.122p . (|4.113p . (|4.104p and (j4.123p we deduce that for n < n* 

\\s"''\\<Cxi\\e"\\ + m), 

and also 

||e"''||i<CA. 
We consider now e"'^. Since 

^",2 _ ^n,2 ^ ^„,2 ^ |(W.",2^,,2 _ [/-2^n,2) ^ 1 (v^n,2^„,2 „ ij;^.^^;^.^) 
= (M/«^2gn,2)^ _ gn,2gn,2 ^ 



(4.123) 
(4.124) 

(4.125) 
(4.126) 



and 

we have by (|4.117p for n < n* 



T/n,2T^n,2 TTn,2 TTn,2 (-\rn,2^n.2\ ^n,2^n,2 
^ ^ ^hx ^ ' £ ' )x- £^ e ' , 



1F"^'-f"^'||-i<Ca(|1£"'2| 



and from (|4.116l) and (|4.119p . for n < 



(4.127) 



(4.128) 



||F"''-K''||-i<CA(|k"|| + ||e"||i). 
Therefore, we conclude by ()4.100p . for n < n* 

\\e-^%<C\{\\e-\\ + \n). 

We have now reached the final stage of the estimation of the errors of the fully discrete scheme. First, we 
estimate e"+^: Since 

pn,3 _ pn,3 ^ ^„,3 ^ 3 (yf^«,3^n,3 „ U^^'^'Ul^) + ^iV^'^V'^^'^ - Kx) > 



and 



T/n,3Trn,3 rjn.^ Tjn,3 /T/n,3_n.3\ _n,3_n,3 
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we have from (I4.126P for n < n*, 

\\F-''-F-^'\\^,<C,i\\e--'\\ + \\e-'%), 
and therefore, from (|4.125l) . (I4.128P 

\\F--^-F-'^^,<Cx{\\s-\\ + \\e-h). 
Thus, by (|4.102p . (|4.109p . (|4.118p . (I4.127p . and above inequahty, we obtain for n < n* 

||e"+i||i < ||e"||i + C.MIk"ll + lle"||i) + , (4.129) 

which is the required recursive inequahty for ||e"+^||i. To obtain an analogous relation for requires 



more work. Observe that 
and 



n,3 _ ff,n,3 _ ^n,3 j_ 1 (-yn,3^n,3 _ ^n,3jjn,3^ 



By (|480)) and (I4.124P we obtain 



Thus 
where 



(Rn,3 rf)"'3 Ijjn n k j jn n , 02^7 jn n aia2k^ TTn„n , ~n 

^ — 2^ ^ ~ T,^ Pi ^ P2 16 ^ P3+74) 



Hence, by (|4.104p . (|4.113p . and (|4.132p we obtain 



where 



74" = ^(74"). , P4 = P{U''Pl). ■ 

By (|4.130p and (|4.126p we obtain for n < n* 

Il74"lli < CA||e"'3|li + Cfc|l73"lli + CAfc||e-3|li , 
from which, using (I4.128p . (|4.123p and (|4.125p we see that for n < n* 

\\l2\\i<Cxi\\e-\\ + n), 

and 

l!7?ll <C^A(||£"|| + ||e"||i). 

From (|4.103l) and (|4.112p we obtain 

- '^>l), + 62P($"'i - <Pl% = 1(^1 + fe2)pr - + fci7i + ^^272 

= 4P1 - 24P2 + 67l + 372 ■ 

From (I4.121P and (|4.13ip we see that 



(4.130) 

(4.131) 
(4.132) 



(4.133) 



-^nr^/'4+&373"+fe474" 
iPl - V2P2 + 48/^3 - 2406/^4 + 373 + 674 ■ 



Therefore, in view of ()4.10ip . we have 



where 



75" = b(7i"+74") + |(72"+73"), 



(4.134) 



satisfies, in view of (14.1061) . (14.1141) . (|4.123p and (|4.133p . tfie inequality 

||7?ll<CA(||£"|| + ||e"||i), (4.135) 
for n < n* . We begin now the estimation of 

^" := - fp? + f P2 - + 2rmP4 ■ (4-136) 

We have 

= - fc(£^ pl^) + ^ [llprf + (e^ p'^')] - 1^ [(£", P3 ) + 3(P^ P^)] 

„n 1 1 2 



^[ll/'2lP + + |(Pi,P3)] - me [(Pl.P?) + 2(P^,P?)] 



' 64 

+ ^[llPalP + UP2,Pi)] - 2#16(P3\PJ) + OTil/'Jr ' 

or 

||^n|l2 \\ n\\2 L -f " i fe^ fe'' f'l i fe'' fn i fc'^ J"ri k' fn , fc*^ ll^"l|2 (A 1 

llf'^ II — F II ^ '^/l + —J2 ~ 24-'3 + 64''4 ~ lAT&h + 48^/6 ~ 242-16 J 7 + 24^06^11/5411 ' (,4.i^3^J 

where the numbers /",...,/" have been defined in the obvious manner. By (|4.104l) we have 

A" = (e'NpJ') = -(C/"4\£") = i(C/re",£") , 

and therefore 

l/ri<C||e"|p. (4.138) 

By (|iTTO)) and 

(e", p^) = -(C/"4\ p^') = -((C/"e")., p?) + (f/^e", = -||p'i 11' + (t^^^", P?) • 

Since 

/2' = I|P?II' + (£",P2), 

we conclude that 

/2" = (C/,"£",p?)<f||£"|p. (4.139) 

Now 

In addition, using (|4.122p . (|4.104p . and (|4.113p we have 

(£",p^') - -(c/"4\P2) = ~{pi,P2) + (t^r£",P2) = -Kf/.>r,p?) + me\p-^) . 

Hence, since 

/3" = (e",P3) + 3(pr,p^), 

we have 

73 = (C^:>; Pi > Pi ) + (C^x e > P2 ) , 

from which 

\m<mn'- (4.140) 

Now 

(p?,p^) = -(C/>L,P2) = -((f/>?).,P2) + (0?,P2) = -IIP2IP + (0?,P2), 

and by (|4.132p 

(£",P4") = -(C/"£",P3) - -((C/"£").,P?) + (C/."£",P3) = -(pr,P3) + ([/."£", P3) 
= IIP2f -(t^rP?,P2) + (C/:e",P3)- 

So, since 

/4" = l|P2lP + i(e'\P4) + |(Pl,P3), 

we have 



h = [U^Pi,P2) + 3(U^e ,^3), 

and thus 

l/ri<fT||£"f . (4.141) 

We also have 

(P2,P3) = 5(t^"P2,P2), 
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and 

Since 
we have 

/5" = i(02,P2) + (Ol,P3), 

and therefore 

\m<§.\\e'T. (4.142) 

Similarly, 

/e =IIP3ir + i(P2,PJ), 
/e" = i(C/>2,P3) - hMf < - ^llPsf • (4.143) 

l/"l<|^lk"f ■ (4.144) 



Since 
we have 
Finally, 
and thus 



/l6 

Now, from the inequalities (|4.138l) - (|4.144|) and the fact that 



\p7\\ < f IIP3' 



where C is a constant independent of k and h, we obtain in (I4.137P 

< (1 + c,k)\\s^r + 2T^(# - i)iiP3f . 

Thus, for A < Aq = 32/C we will have the following estimate of ||cr"|| 

lk"||<(l + C.fc)||e"||. 

Therefore, from (|4.134p . (|4.136p and (|4.135p . we finally obtain 

lk"+^|| <||£"||+C.fc(||£"|| + ||e"|li) + 115^11, 
for n = 0, 1, . . . , n*. From this inequality, (|4.129p and Lemma 4.5 we obtain for n < 



n+l 



II + ||e"+^||i < (1 + CAfc)(||£"|| + ||e"||i) + Ckik^ + h''-'^y/ha/h) . 



Therefore, by Gronwall's Lemma, and (|4.69p we see that 

||e"|| + ||e"||i <CA(fc* + /i3-^v/ha77^), 

for n — 0,1, . . . ,n* . Hence, if k and h are sufficiently small, the maximality of n* contradicted and we may 
take n* = M — 1. The conclusion of the proposition follows. □ 

Remark 4.5. Arguing as in Remark 4.1 we have again ||C/," — u{t"-)\\oc = 0{k^ + h'^-^^/ln l/h). 

Remark 4.6. The analogous results of Lemma 4.5 of Proposition 4.3 hold for the system (jCBp as well. 

5. A NON-STANDARD GALERKIN METHOD 

In this section we analyze a semidiscrete Galerkin approximation of (jGBp and (|SGBp which is non-standard 
in the sense that it approximates ry by piecewise linear, continuous functions and u by piecewise quadratic, 
Cq functions, defined on the same mesh. (The method may be generalized using the analogous pairs of 
higher-order splines but here we confine ourselves to the low-order case.) In the analysis, we consider only 
the uniform mesh case, where a series of superaccuracy results for the error of the L^-projection onto the 
space of piecewise linear continuous functions affords proving optimal-order error estimates. 
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5.1. Preliminaries and the basic results. Using the notation of paragraph 2.2, in addition to Sf^ we let 

{^^C'- 0|[,^,,^.^,] e P2 , 1 < J < A^, 0(0) = 0(1) = 0} . 

In this section we shall denote by a(-, •) the usual bilinear form a(^/;,x) = (0, x) + §(0''jX') defined for -0, 
X G S'j^ ^^"^ • ^0 ~^ elliptic projection operator defined for w € -ffg by {{Rhv)', x') = 

(w',x') Vx € S'^^Q. The system 0CB^ is discretized as follows: Seek rjh : [0,T] ^ S^, uu : [0,T] ^ 5*^ g such 
that for < t <'r 

{Vht,(t>) + [uhxA) + \{{VhUh)x,<f)) , yi'^Sl (5.1) 
a("/it,x) + (»7fti:,x) + f (w'i'«/i2;>x) + 5(%?7/ix,x) = 0, Vx e 5'^,o (5-2) 

and 

f]h{0)^Prio, Uh{0)^RhUo, (5.3) 
where P : S^^ is the L^-projection operator onto S^. (We will treat only the case of (|SCBI) . that of 

(jCBp being similar.) The o.d.e. initial-value problem represented by (|5.ip - (|5.3p has a unique solution for any 
T < 00. In fact, taking = rjh in (|5.ip and x = in (|5.2I) and adding the resulting equations, we obtain 
for t > 

\Mt)f + \\u,M\l = \MO)f + \\uMl, (5.4) 
and the desired conclusion follows by standard o.d.e. theory. 

The plan of the section is as follows: We will start by stating the preliminary Lemmas 5.1-5.4 with the 
aid of which we will prove our main error estimate. Theorem 5.1. The proofs of the Lemmas depend on 
some superaccuracy estimates of the projection error v — Pv for smooth functions v; these are proved in 
§5.2. (Lemmas 5.5-5.10). Recall first the following standard estimates, [19j,[3(i,. Suppose that 77 S and 
u G Cg. Then, if p = 77 — Prj, a = u — RhU, we have 

WpW+hMiKCh" , \\a\\ + h\\a\\i < Ch\ ||a||oo + /i||p||oo < C/i^ 

Lemma 5.1. Let rj G v E Cq. If p — t] ~ Prj and is such that 

then IICII < Ch^. 



Proof. The result follows from Lemma 5.9 and Lemma 2.1(m). □ 
Lemma 5.2. Let rj e and u e Cq . If p ~ rj ~ Prj and ( £ is such that 

then IICII < Ch^. 

Proof. The result follows from Lemma 5.10 and Lemma 2.1(m). □ 
Lemma 5.3. Let rj E and u E Cq. If p — rj — Prj, a ^ u — RhU, then 

(i) (p',V) = y^jesi„ 



(ii) (a', 0)^0 V0e5i 



Proof. The identity (z) follows from the observation that (p', -0) = (p, -0') = 0, since £ Sf^Q implies ip' e 5*^. 
In order to prove (m), we use an argument similar to that of Wahlbin, [30, p. 5]. Let be the usual hat 

function basis of S'^ associated with the uniform mesh Xi = (i — l)h, 1 < i < N + 1, i.e. with the property 
that (jjiixj) = Sij. Then, the function ipi defined for 1 < i < iV + 1 on [0, 1] by "01(3;) = Jq <Pid£, — x (j)idS_ 
belongs to and satisfies — (jji — Jq (t^id^. Since (cr', 1) = 0, it follows that {a',(l)i) = [a' — 0. 
Hence, (ti', 0) = for any cj) G and the proof of (ii) is completed. (Since Pu' = (Rhu)' as it may be easily 
established, (ii) implies that cr' = it' — Pu' .) □ 

Lemma 5.4. Let rj G and u £ Cq. If a = u — RhU and ( € is such that 

(C,0) = ((w)',0) V0e52, 

then IICII < Ch^. 
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Proof. Consider 7^ = ((rjcr)' , Then, by the final remark hi the proof of Lemma 5.3, 

7, = (7/cr, 00 + (r/cr', = (?7'fT, 4>,) + (77(m' - Pm')>.) ■ 

Now, since |1(t|1oo < Ch^, we have that \{7fa,(j)i)\ < C\\a\\ 00 Jq (t'idx < C/i''. In addition, by Lemma 5.10 
we obtain that \{tj{u' - Pu'),(t>i)\ < Ch^. We conclude that I7I = (YlfJi^ ifY^^ < C'^^■^ from which the 
conclusion of the Lemma follows in view of Lemma 2.1(m). □ 

With the aid of these lemmas we may establish the basic result of this section: 

Theorem 5.1. Suppose that the solution {ri,u) of ^SCB\ is such that 77 e C(0, T; C*), r/t G C(0,T;C^), u, 
Ut G C(0,T;Cq) and let {rjh,Uh) be the solution of the semidiscrete problem 1^5.1^ - ^5.3^ . Then 



(i) gmax^dl^W - + f'Mt) - Vh{t)\\) < Ch 



3 

3 



(m) ^max^(||u(t) - u/i(t)||oo + hWvit) - Vh{t)\\oo) < Ch 

Proof. Let p ~ rj — Prj, 6 = Prj — rjh, a — u — RhU, and ^ ~ RhU — Uh. Then rj — t^h — p + OjU — Uh = o' + ^ 
and, by Lemma 5.3 

i9t,(f>) + {^^,^) + ^{{vu~VhUh),,^)^0 yep ^ Si, (5.5) 
a(6, V") + (^x, V') + - UhUhx,^) + ^{vVx - rihVhx,il^) = ~i<yt,ip) VV- G S^ Q , (5.6) 

for i e [0,T]. Since 

T]u - r7hM/i = r]{u - u,,) + ^(77 - r]h) - (77 - 77,i)(u - Uh) = rj{(7 + £) + u{p + 0) - {p + 0){a + ^) , 
it follows that 

■qu - TjhUh ^ria + up-9^ + f , 
where f — rj^ + u9 — pa — pS^ — da. In addition, 

VVx - VhVhx = viVx -Vhx) + Vxiv - rih) - {v - 'nh){rix - rjhx) 

= vipx + Ox) + vx{p + o)-{p + e){px + Ox) 

= {ilP)x + {Tld)x - {pO)x - PPx - d9x , 

whence 

vvx ~ mmx = {'np)x - oo^ + gx , 

where g — rjO ~~ p6 ~ ^p^. Finally, 

UUx - UhUhx = u{Ux - Uhx) + Ux{u - Uh) - {u - Uh){Ux - Uhx) 

= u{ax + ^x) + Ux{a + i) - {a + i){ax + ix) 
= {ua)x + {uCjx - {aC)x - aax - £,£,x , 



I.e. 

where f — ua + — a£^ — i(T^ — i^^. Thus we rewrite (j5.5p and (j5.6p as 
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UUx - UhUhx = fx 

we rewrite (j5.5p £ 

0) + i((w)., + UMx,4') - hm)x,^) + i(/.» = V0 G ^^ 



Putting = 0, X = ^ in the above and adding the resulting equations, we have for < t < T 

Hm\' + iieii?) + Hvcr)x,9) + ^iiup)x,9) + ^ifx,9) 

+ ^{{vp)x,0 + lifx,0 + H9x,0 - • 

We estimate the various terms in (|5.7p as follows. By Lemma 5.4 we have 

\{{w)x,o)\<ch^\\e\\, 

and by Lemma 5.1 

\{{up)x,e)\<ch^\\e\\. 
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(5.7) 



Using integration by parts and the standard estimates for p and cr we also have 

\{f.„e)\<cu\\i\\0\\+ch^\\0\\+c\\e\\\ 

In addition, from Lemma 5.2, since e Sf^ 

\{{rip)..0\ = \{riP.£..)\<Ch^m\i- 
Finahy, it fohows from integration by parts and standard estimates for a that 

\{L^)\^\{li.)\<ch'm\i + cm\\ 

and by standard estimates for p that 

\{9.,o\ = \{9.^.)\<c\\emh + ch'm\,. 

Hence, from (|5.7p and these estimates we deduce for < t < T 

~m? + m\\)<c{h' + \\er + m\\), 

and, consequently, by Gronwall's Lemma and (j5.3p that 

\\e\\ + U\\i<ch^. o<t<T, 

and the conclusions of the Theorem follow in view of the approximation properties of S'^ and Sf^ q, Sobolev's 
theorem and the L°°-L^ inverse inequality in Sf^. □ 

Remark 5.1. In view of Theorems 1.3.1 and 1.3.2 of the conclusions of Theorem 5.1 remain valid if in 
place of RhUo we use the elliptic projection of uq that is induced by the bilinear form a(-, •) as was done in 
previous sections. 

5.2. Superaccuracy estimates for the error of the projection. The proofs of Lemmas 5.1, 5.2, and 
5.4 depend on some superaccuracy properties of the error of the projection onto Sf^ of a smooth function, 
in the case of uniform mesh. 

Suppose that rj £ , u £ Cq , and let P be the projection onto S'^ and p = i] — Prj. To estimate 
(wp, 0^), where {(pi] is the standard basis of S*^, it suffices to estimate the integrals 



updx , 

where Ii = {xi, Xi^i), i ~ 1,2, . . . , N . Putting 2^^+1/2 = {xi + .Xi+i)/2, wc will have 

updx = u{x,+i/2) / pdx + u'{xi+i/2) / {x - Xi+i/2)pix)dx + ^ / {x~'Xi+i/2fpix)u"{yi)dx, 



for some yi — yi{x) G li. Since ||p||oo = 0{h'^), the third integral in the right-hand side is of 0{h^). Now, if 
X £ h 

P{x) = p{Xi+l/2) + {x- X,+i/2)p'{Xi+i/2) + ^{X- X,+i/2fv'\x,+i/2) + ^(x - X,+ i/2fv"' in) , 

for some r,; = Ti{x) G /,;. Hence 

y (x - x,+i/2)p{x)dx = p'(x,+i/2) y (2; - x,+i/2fdx + 0{h^) = j^p' {x,_^_i/2) + 0{h^) . 

Thus 

j updx = u{x^+i/2) j pdx + ^p'{x,,+i/2)u'{x.,+i/2) + 0{h^) , 

and, consequently, we must estimate the integral of p on /; and the derivative of p at the midpoint of 
Both of these quantities turn out to be superaccurate as we prove by elementary techniques in the following 
two lemmas: 



Lemma 5.5. Let -q e . Th 



max |p'(a;,+i/2)| < Ch^ 

l<i<N 
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Proof. For x € li we have 



p'ix) = r^'ix) - {PvYix) = rj'ix) - ^1±1^ 



I.e. 



p'(.X^+l/2) = v'{x,+ ^/2) - ^^±lj-^ , 1<1<N, (5.8) 

where Ci is the ith component of the solution c of the hnear system Gc — b, where Gij = (jfi), 1 < i,j < 
N + 1 and bi = (77, i = 1,2, . . . , N + 1. The equations of the system may be written exphcitly as 

2ci + C2 = 661//1 , 

ci + 4c2 + C3 = 6b2/h, 

Ci-i+4ci + Ci+i=6b,/h, i^3,i,...,N -1, (5.9) 

CAr_i + 4cAr + Cn+1 = Qb^/h , 

cn + 2cAr+i = 6bN+i/h . 

Hence 

3(c2 - ci) + (C3 - C2) - 6(62 - 26i)//i , 
(ci - Cj_i) + 4(cj+i - Ci) + (cj+2 - Cj+i) = 6(6j+i - 5i)//i, 2 < i < iV - 1 

and 

(cat - CN-i) + HcN+1 - Cn) = 6(2feAr+i - bN)/h. 
Therefore we obtain the following equations for the differences (c^+i — Ci)/h: 

3(c2 - c,)/h + (C3 - C2)/h = 6(62 - 26i)//i2 , 
{ci - Q_i)//i + 4(q+i - Ci)/h+ (cj+2 - Ci+i)/h = 6(6i+i - 6i)//i^ , i ^ 2,3, . . . , N - 1 , 
{cn - CN-i)/h + 3{cN+i - CN)/h = 6{2bN+i - bN)/h^ ■ 

Hence, if we put e[ :— ^'(2^^+1/2), i — 1, 2, . . . , A^, then (|5.8p implies that e' = (e^) is the solution of the linear 
system 

Ae'^r', (5.10) 
where A is the N x N tridiagonal matrix with An = Ann = 3, An =4, 2<i<iV — 1, and Aij = 1 if 
|z — j| = 1, and r' — {r[,r2, ■ ■ ■ , r'pf)'^ , where 

r[ = 3r?'(x3/2) + v'{^5/2) - 6(62 - 2bi)/h^ , 

r- = r;'(x,_i/2) + 4r7'(a:,+i/2) + - 6(6.+i - b,)/h^ , i = 2,3, . . . , N - 1 , (5.11) 

i^'n = v'{xn-i/2) + 'W{xn+i/2) - 6(267v+i - bN)/h'^ ■ 
We will now show that there exists a constant C depending only on ||?7"'||oo, such that 

max Ir'l < Ch^ . (5.12) 

l<i<Af 

To estimate r[ note that 

62 — 261 = / rj(f)2dx — 21 i](j>idx = ^ h)ridx — ji {x — h)rjd: 

J huh Jh Jh Jin 



Tjdx 

/1U/2 J Ii J Ii J I2 -I huh 

^- j;Ji - tJ2 + J3 ■ 



Using Taylor's theorem we have 

Ji= [ {x-h)r]dx = -!^rj{h) + !^T]'{h)-i^Tj"{h) + 0{h^). 
Jh 

J2= / (x-h)r]dx = i!^rjih) + !^ri'ih) + l^rfih) + 0{h^), 



J3= / 7]dx^2h7]{h) + ^7]"{h)+Oih^), 



huh 
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and thus 
In addition, 



Therefore, by (|5.1ip . we see that r'l — 0{h^). For rjy we have 

26jv+i — bjy = 2 / ri(f)N+id,x — / ri(f>Ndx 

J In ^ In -i^In 

[x — X M)Tidx — I [x — XN)'ndx - 



In- 



•qdx 



In-i^In 



By Taylor's theorem 



J 2 = 



3 



,h = 2hr^{xN) + ^v"{xn) + 0{h^) 



and thus 
Hence,by ((5lT1) 
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7V+1 — Ojv 



= (2:^-1/2) + 'i-n'{xM+i/2) - ^{'2-bN+i - bN) 
= iv'ixN) + hri"{xN) - ^v'{xn) - hri"{xN) + 0{h^) = 0{h^) . 
For r^,2<z<A^ — 1, we have 



bi+i -bi ^ r](l)i+idx - / r]<j)idx 

'li'Jli + l Jli-iUh 



h / {x - x,+i/2)iida 



[x - Xi+3/2)r]dx - \ I (a; - x,_i/2)?7d; 



+ 



■/ ridx~\ 2j^_ IJ3+ IJ5. 

Jii+i Jii-i 



By Taylor's theorem we obtain 

Ji = ??'(.x,+i/2) J^ix- x,+i/2?dx + 0(/i5) = f^?7'(a;.+i/2) + 0(/i^) . 

Similarly 

J4 = /ir/(a;,+3/2) + |^?7"(a:»+3/2) + 0{h^) , 
J5 = h'q{x,_i/2) + ^v"{x,^^/2) + Oih") . 

Since now 

J4 ~ J5 = 2/i2?7'(x,+i/2) + 0{h^) . 
Substituting (|5.14p - (|5.17p into (|5.13p we obtain 

bi+i -b, = ^v'ixi+1/2) - ^'n'ix,+3/2) - ^v'ixi-1/2) + 0{h^) , 

and finally, by ([ETT|) 

r', = ?7'(a;j-i/2) + 4?7'(a;j+i/2) + 77'(2;-*+3/2) - -^{bz+i - 6i) 

= i('7'(a:»+3/2) - 2v'ix^+l/2) + v'ix^-l/2)) + 0{h^) = 0{h^) . 



(5.13) 



(5.14) 

(5.15) 
(5.16) 



(5.17) 



We conclude therefore that (|5.12p holds with a constant C — C(|l77"'|loo)- Writing the linear system (|5.9p as 
\Ae' = \r' with \A = I -B, where B is a N x N matrix with ||B||oc = 1/2. Therefore, ||(/- B)-i||oo < 2 
and consequently 

max le'l < i max Ir'l < Ch'^ . 

l<i<N I - l<i<Ar I «i - 

□ 

Lemma 5.6. Let ?/ G C"^. Then for some constant C = C'(||'?"'||oo) we have 

max I / pdx\ < Ch"^ . 

l<t<N^Jj^ I 

Proof. We have 

J pdx = y" rjdx - |(cj+i + Ci) , i = 1, 2, . . . , , 

where Ci is the i-th component of the solution c of the linear system Gc — b with G and b as in the previous 
Lemma. From the equation (|5.9p we see that 



5 • |(ci + C2) + |(C2 + C3) = 3(25i + 62) , 
|(c,_i +c,) + 4- |(c, + c,+i) + |(c,+i +c,+2) = 3(5,+i +6,), 2<i<N-l, 
|(cAr_i + Civ) + 5 • Kctv + cjv+i) = 3(6jv + 26w+i) . 

Hence, if 



Ei := J pdx , i = 1, 2, . . . , , 

then e = (ei, £2, . . . , ew)^ is the solution of the system 

re = r, (5.18) 

where T is the tridiagonal N x N matrix with Fn = Fjvat = 5, Ta — 4, i — 2,3, . . . , N — 1, and F^ = 1, 
when |i — j| = 1, and where r = (ri, r2, . . . , r^)'^ with 

/2 



li- 



ri = 5 / ?7(i.T + / i]dx — 3(26i + 62) , 
J12 

T]dx + -i ridx+ ridx - 3{bi+i + bi) , i ^ 2,3, . . . , N - 1 , (5.19) 
rjsi ^ I rjdx + 5 rjdx — 3{bN + 2bN+i) ■ 

J In -I J In 

We will show that = 0{h'^), for i = 1, 2, . . . , iV . 
For ri, note that 

261+62 = 2 / rjipidx + / ri(f>2dx = / rjdx — / {x — h)rjdx . 
Jii J11UI2 J11UI2 JI1UI2 

Hence, by Taylor's theorem 

26i + 62 = 2hrj{h) - ^r]'{h) + i^f]"{h) + 0{h^) . (5.20) 

Also 

rjdx = hr^ih) - l^v'ih) + + 0{h^) , 



T]dx = h7]{h) + ^T]'ih) + \-T]"{h) + 0(/i4) . 

I2 

From these relations and (15.191) and (|5.20p we obtain 

n = b{hij{h) ~ i^rj'ih) + !^rj"{h)) + (hrjih) + {h) + '^^\h)) 
- 3{2hrj{h) - ^v'{h) + !^r]"{h)) + 0{h^) = 0{h^) . 
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For we observe 



^AT + 26jv+i = / ri(f>Ndx + 2 / -qcfyisi+idx = / [x — x isi)i'ldx + I rjdx . 



Hence, by Taylor's theorem 



bN + 26^+1 = 2/i77(a:Ar) + ^ri'{xN) + ^r;"(a;Ar) + Oih") . (5.21) 
Also 

■qdx = hr]{xN) + ^rf{xN) + ^r]"{xN) + 0{h'^) . 

From these relations and (|5.19p and (|5.21l) wc sec that 

rN = {hri{xN) - ^t]'{xn) + ^ri"{xN)) + h(hrj{xN) + ^r]'{xN) + ^v'ixN)) 
~ 3{2hTj{xN) + ^v'{xn) + 'tv"{xn)) + 0{h^) = 0{h^) . 
For Ti, i = 2, 3, . . . , — 1, we have 



bi + bi+i = / ri<j)idx + / rjc/jidx 

= y {x ~ Xi^i)ridx + J rjdx + j {xi+2 - x)r]dx . 



+1 



Thus, by ^^IM wc have, for 2 < i < A - 1 



r]dx~f {x - Xi^i)r]dx - f {xi+2 - x)r]dx =: Ji - fJ2 ~ j^Js ■ (5.22) 

/i_iU/iU/i+i -Jli-i Jli+1 

By Taylor's theorem 

J2 = f r,(a;,_i) + f r7'(^,-i) + ^^"(a^.-i) + 0{h'>) , 

i.e. 

J2 + h = f (?7(^»-i) + v{x^+2)) + f (r?'(:E.-i) - v'{x,+2)) + ^(?7"(^»-i) + r;"(a;.+2)) + 0(/i') 

= /i'»7(a^«+i/2) + ^-n^x^+m) + o(/i5) . 

In addition, 

Ji = 3/iry(x,+i/2) + \^Tl"{x,+y2) + 0(/i*) . 
Substituting in ([Q^]) and in (|5?T9l) we see that for 2 < i < A^ — 1 

n = 3hrjix,+,^2) + ^v"{x,+,/2) - 3h7j{x,+,/2) - 2|lr;"(x,+i/2) + 0{h*) = 0{h*) . 
Hence, there exists a constant C depending only on ||77"'||oo such that 

max \r,\ < Ch^ . (5.23) 

From (|5.18l) we obtain that jTe = |r, where jT — I ~ E, and i? is a A^ x A^ matrix with ||-E||oo = 1/2. 
Hence - E)~^\\^ < 2 and thus 



max leil < i max IrA < Ch^ . 

l<i<N 1*1-2 i<i<^ I 'I - 



□ 



If we now assume that rj is smoother, then, using a result by Demko, [T7\, (see also de Boor,, 16 ), we 
may show that Jj pdx ~ 0{h^) provided that the endpoints x^, Xij^i of the interval Ii are at a distance of 
0{h\x\ 1/h) from the boundary of [0, 1]. 
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Lemma 5.7. Let rj G C*. Then, there exists a constant C such that if 

dist{xi,0) > Chlnl/h and dist{xi^i,l) > Chlnl/h , 

then 

pdx = 0{h^) . 



Proof. We will first show that due to our hypothesis of increased smoothness of 77, we now have 

r, = Oih^), i = 2,3,...,N-l, 

where the were defined in (j5.19p and written in the form (|5.22l) . Since 77 G C^, using Simpson's rule we 
approximate the integrals in (|5.22p by 

Tjdx = ^ + Mx^+l/2) + V{x^+2)) + , 

/i_iU7iU7i + i 



X-—Xi 

h-1 

^^Tjdx = |(r;(a.,+i) + 2r7(a;,+3/2)) + 0{h^) , 



+1 



yielding 

n = ^{v{xi-i) - 2ri{x,_i/2) - + 477(x,+i/2) + ~ '^v{xi+3/2) - v(xt+i)) + 0{h^) . 

Hence, using Taylor's theorem 

n = |(^r7"(x,_i/2) - 'i^i'{x^+l/2) + ^77"(a;,+3/2)) + 0{h^) 
= l^{n"{x^,_„2) - 2r,"(a;,+i/2) + v"{x^+m)) + Oih") = Oih^) . 
We conclude that for constants Ci, C2 independent of h we have, in view of (I5.23p . that 

In I < Cl/^^ max \r,\ < Cs/i^ , |rA,| < Cih'^ . (5.24) 

2<i<Af-l 

We also recall from Lemma 5.6, cf. (I5.18p . that the Si :— Jj pdx, I < i < N, form the solution of the system 

lTe={I-E)e=lr, 

where F was defined after (|5.18p . Hence 

N N-1 

= Y.i'^-'hr, = (r-i).iri + ^ (T-%r, + (T-^^rN . (5.25) 

Now, using Proposition 2.1 of [T7], we see that there is a constant C3 such that 

|(r-i),i| <C32-l^-i|, and |(r-i).Ar| <C32-I^-^I. (5.26) 
Therefore, if dist{xi,0) > -y^hlnl/h and dist{xi^i, 1) > j^ft-lnl//i, there follows that 

i-l>j^lnl//i and iV - j > In , 

i.e. 

2-l^-i|</i and 2-l'-^l</i, 

and by ((5:26| 

|(r-i).i| <C3/i and \{V-%N\<C3h. 

Hence, for the indices i for which 

l + _i_lnl//i<i <7V- i^lnl//!, 
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we will have from (|5.25l) and (I5.24p that 



Af-l 



i=2 

AT 

< CaCi/i^ + max V \(^-^)^AC2h'^ + CgCi/i^ < C/i^ 

1 <i<-N ^ — ' 



l<i<N ■ 



□ 



We now prove, as consequences of these results, several superaccurate estimates on weighted integrals of 
the error of the projection. 

Lemma 5.8. Let rj E C*, u E Cq, and p = -q — Prj. Then, 

max I / updx\ < Ch^ Inl/h . 

Proof. If 1 < i < iV we have for ti = ti{x) E li 

Updx ^ u{Xi+i/2) I pdx + u'{Xi+i/2) / {x ~ Xi+i/2)pdx + ^ / {x~Xi+i/2Ypu"{ti)dx. 



The third integral of the right-hand side is of 0{h^). For the second integral, by Taylor's theorem we have 

^ (x - x^+ii2)pdx = Y^p'(2:^+i/2) + 0{h^) , 

and using Lemma 5.5, 



We conclude that 



If now i is such that 

dist{xi^Q) > Ch\x\l/h , and dist{xi+i,l) > Ch\\il/h , 
where C is the constant in the statement of Lemma 5.7, then 

updx^O{h^). (5.28) 



(a; — Xi^i/2)pdx = 0{h^) . 
updx — u{xij^ii2) I pdx + 0{h^) . (5-27) 



If dist{xi,0) < Chlnl/h, we have 

""(2^1+1/2) = u{xi) + !^u'{n) , u{xi) = Xiu'{yi) , 
for some E (xi, 2:^+1/2) and yi E (0,2;^). Therefore 

1^(2:2+1/2)1 < Ch\nl/h, 
and from Lemma 5.6 and (j5.27p we have for such i that 



updx\ < Ch^lnl/h. 



(5.29) 



If finally dist{xi-i-i, 1) < Chlnl/h, we have 

u{X'i.+i/2) = -"(2^1+1) - ^u'{t,) , u{xi+i) = (1 - Xi+i)u'{yi) , 
for some Ti E (2:^+1/21 2:^+1) and yt E {xi+i, 1). Hence 

l"(a;i+i/2)| < Chlnl/h, 
and it follows again from Lemma 5.6 and (|5.27p that in this case too 

1^ updx\ < Ch^lnl/h. 



(5.30) 
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The conclusion of the Lemma fohows from (|5.28p - (|5.30p . □ 
Lemma 5.9. Let rj G C*. u E Cq, p = i] — Pi], and 

I3,:={up,(f>'^, i = l,2,...,N + l. 

///3 = (/3i,^2,...,/3Ar+if then 

\P\<Ch'-', 

where \ ■ \ is the I2 norm in MJ^^^ . 
Proof. We have 



/3i = -i / updx , 
Jh 

Pi — jl j updx — \ j updx , i = 2, 3, . . . , A'^ , 



Pn+1 = i updx . 

J In 



If dist{xi^i,Q) > Chlnl/h and dist{xi+i,l) > Chlnl/h, i.e. ii i E K, where K are the integers in the 
interval [2 + C In - (7 In , then by ((OS)) 

m<Ch\ (5.31) 

If on the other hand dist{xi^i,0) < Chhil/h or dist{xi+i,l) < Ch\nl/h, i.e. if i G Ki, where Ki are the 
integers in [1, 1 + Cln U [N + 1 - C \nl/h, N + 1], then, by (lOO)) and ((OOl 

lAI < Ch'^lnl/h. 

Therefore, 

N+l 

< {I + 2C\nl/h)Ch^{\nl/hf + C{N - 1 - 2C\nl/h)h^ < Ch'' . 

□ 

Lemma 5.10. Let w G C^, w G and 

P, = {v{w-Pw),dp^), i = l,2,...,iV + l. 
///3 = (/3i,/32,...,/3iv+i)^, then 

|/3| < C/z3-5 , 

where \ ■ \ is the I2 norm in M.^^^. 
Proof. For a; G Ii we have 

v{x) = v{h) + 0{h) , 

and since — Pw||txj = 0{h?), 

v(x){w - Pw){x) = v{h){w - Pw){x) + 0{h^) , 

and therefore 

Pi = Oih") . 

If a; G /i-i U i = 2, 3, . . . , iV, we have 

v{x) — v{xi) + 0{h) , 

whence 

v{x){w - Pw){x) = u(a;i)(w - Pw){x) + OQi^) , 

and so 

P, = 0(/i^) . 

Similarly, Pn+i = 0{h^), and finally \P\ = 0{h^-^). □ 
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5.3. Numerical experiments. We considered first the (ICBI) system, that we discretized by the nonstandard 
method analyzed in the previous two sections using the subspaces Sf^ and S'j^ q for approximating 77 and u, 
respectively. (We considered the nonhomogeneous system with a suitable right-hand side so that the solution 
is 77 = exp(2f)(cos(7ra;) + .x^ + 2), u — exp(a;i)(sin(7ra::) ^ — x^).) We integrated the semidiscrete system 
up to T = 1 by the fourth-order accurate explicit Runge-Kutta method analyzed in paragraph 4.3 with time 
step k — /i/lO; the temporal discretization error was negligible in comparison to the spatial error. Table 



N 




order 


u 


order 


40 


0.1250(-2) 




0.4057(-5) 




60 


0.5555(-3) 


2.001 


0.1199(-5) 


3.008 


80 


0.3124(-3) 


2.000 


0.5051(-6) 


3.004 


100 


0.1999(-3) 


2.000 


0.2585(-6) 


3.002 


120 


0.1388(-3) 


2.000 


0.1495(-6) 


3.001 


140 


0.1020(-3) 


2.000 


0.9416(-7) 


3.001 


160 


0.7810(-4) 


2.000 


0.6307(-7) 


3.001 


180 


0.6170(-4) 


2.000 


0.4429(-7) 


3.001 



N 


n 


order 


u 


order 


40 


0.3342(-2) 




0.5955(-5) 




60 


0.1485(-2) 


2.000 


0.1780(-5) 


2.979 


80 


0.8357(-3) 


1.999 


0.7539(-6) 


2.985 


100 


0.5349(-3) 


2.000 


0.3870(-6) 


2.989 


120 


0.3714(-3) 


2.000 


0.2243(-6) 


2.991 


140 


0.2729(-3) 


2.000 


0.1414(-6) 


2.992 


160 


0.2089(-3) 


2.000 


0.9483(-7) 


2.993 


180 


0.1651(-3) 


2.000 


0.6665(-7) 


2.994 



N 


V 


order 


u 


order 


40 


0.3872 




0.1048(-2) 




60 


0.2581 


1.000 


0.4652(-3) 


2.002 


80 


0.1936 


1.000 


0.2616(-3) 


2.001 


100 


0.1549 


1.000 


0.1674(-3) 


2.001 


120 


0.1290 


1.000 


0.1162(-3) 


2.000 


140 


0.1106 


1.000 


0.8540(-4) 


2.000 


160 


0.9678(-l) 


1.000 


0.6538(-4) 


2.000 


180 


0.8603{-l) 


1.000 


0.5166{-4) 


2.000 



Table 5.1. Errors and orders of convergence 
cretization, rjh 



e Si, Uh 



e S'jJ Q, uniform mesh. 



CB|) system, nonstandard Galerkin semidis- 



15.11 shows the resulting errors and orders of convergence in various norms. The optimal-order L^- and L°° 
rates of convergence of Theorem 5.1 are confirmed and the rates are as expected. We obtained the same 
rates for the (jSCB|) system. Table 15.21 shows the analogous errors and rates of convergence for the same 



L2. 



N 


n 


order 


u 


order 


40 


0.2214(-4) 




0.1529(-6) 




60 


0.6628(-5) 


2.975 


0.3044(-7) 


3.981 


80 


0.2810(-5) 


2.982 


0.9670(-8) 


3.986 


100 


0.1443(-5) 


2.986 


0.3970(-8) 


3.989 


120 


0.8370(-6) 


2.989 


0.1918(-8) 


3.991 


140 


0.5279(-6) 


2.990 


0.1036(-8) 


3.992 


160 


0.3540(-6) 


2.992 


0.6081(-9) 


3.993 


180 


0.2488(-6) 


2.993 


0.3799(-9) 


3.994 



N 


V 


order 


u 


order 


40 


0.4812(-4) 




0.6323(-6) 




60 


0.1422(-4) 


3.006 


0.1290(-6) 


3.921 


80 


0.5993(-5) 


3.005 


0.4145(-7) 


3.946 


100 


0.3066(-5) 


3.004 


0.1713(-7) 


3.958 


120 


0.1773(~5) 


3.003 


0.8314(-8) 


3.966 


140 


0.1116(-5) 


3.003 


0.4507(-8) 


3.972 


160 


0.7474(-6) 


3.003 


0.2651(-8) 


3.976 


180 


0.5248(~6) 


3.002 


0.1659(-8) 


3.979 



N 


V 


order 


u 


order 


40 


0.5999(-2) 




0.3854(-4) 




60 


0.2654(-2) 


2.011 


0.1152(-4) 


2.978 


80 


0.1490(-2) 


2.008 


0.4884(-5) 


2.984 


100 


0.9520(-3) 


2.006 


0.2509(-5) 


2.988 


120 


0.6605(-3) 


2.005 


0.1454(.-5) 


2.990 


140 


0.4850(-3) 


2.004 


0.9168(-6) 


2.991 


160 


0.3711(^3) 


2.004 


0.6148(-6) 


2.992 


180 


0.2931(-3) 


2.003 


0.4321(-6) 


2.993 



Table 5.2 
cretization, 



Errors and orders of convergence 



Si-, Uh 



g S*^ Q, uniform mesh. 



CBp system, nonstandard Galerkin semidis- 



problem but now discretized so that 77/1 € Sf^, i.e. in the space of quadratic splines, and u/i € S"! g, i.e. in 
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the space of cubic splines that vanish at a; = and at x = 1. We observe e.g. that the associated L^-rates of 
convergence for this higher order pair of subspaces are, as expected, equal to 3 for rj and 4 for u. The (jSCB[) 
system gave the same results. 

The error analysis in paragraphs 5.1 and 5.2 depended strongly on the assumption of uniform mesh. In 
the case of nonuniform mesh we expect some order reduction. For example. Table 15.31 shows the errors 
and orders of convergence that we obtained for (jCBp by approximating (rjh^Ufi) in S'^ x Sf^Q and using for 
both subspaces the exact solution and the nonuniform mesh with which Table [2?3l[a) was produced. Both 
rates are apparently reduced by one now. 



N 


V 


order 


u 


order 


80 


0.1372(-2) 




0.3879(-5) 




160 


0.6849(-3) 


1.002 


0.9569(-6) 


2.019 


240 


0.4564(-3) 


1.001 


0.4241(-6) 


2.007 


320 


0.3422(-3) 


1.000 


0.2383(-6) 


2.004 


400 


0.2738(-3) 


1.000 


0.1524(-6) 


2.002 


480 


0.2281(-3) 


1.000 


0.1058(-6) 


2.002 


560 


0.1955(-3) 


1.000 


0.7775(-7) 


2.001 


640 


0.1711{-3) 


1.000 


0.5952{-7) 


2.001 


720 


0.1509(-3) 


1.000 


0.4702(-7) 


2.001 



Table 5.3. L^-errors and orders of convergence. (|CBI) system, nonstandard Galcrkin 
semidiscretization, rjh € S^, Uh & Sf,^ q, quasiuniform mesh with '^^jl' =1.5 



6. Remarks on standard Galerkin methods for related hyperbolic problems 

We conclude the paper with a section of remarks on the application of standard Galerkin methods on 
some simple first-order hyperbolic problems. Our aim is to show examples of such problems in which the 
techniques previously developed in this paper may be applied to some advantage. 

6.1. Initial-boundary- value problem for a single hyperbolic equation. For T > consider the 
initial-boundary- value problem 

r]t + r)^^0^ 0<a;<l, 0<t<T, 

T]{x,0)^r]oix), 0<x<l, (6.1) 
?7(0,t)=0, 0<t<T. 

It is clear that if 770 € C*^, for integer fc > 0, and tiq ^O) ^ 0, j — 0, . . . ,k, then (j6.ip possesses a unique 
solution 77 given for {x,t) e [0, 1] x [0,r] by 



r]{x,t) 




X > t , 
X <t. 



The solution is classical if fc > 1 and, for example, has the property that r]{-,t) E for all t £ [0,T]. 
For integer N > 2 and h = 1/N we consider the uniform mesh given hy Xi — ih, i — 0,1, ... , N, and put 

Ij = [xj^i,Xj]. We let 

§1 = {0 e C[0, 1] : e Pi, i<j<N, m = 0} , 

and denote by the basis of Sf^ with the property that 4>i(xj) — Sij. The standard Galerkin method 

for (|6.1|) may be then defined as follows: We seek 77^, : [0, T] 5^, such that 

Vh{0) ^Pvo, 

where P is the L^-projection operator onto S'^. Obviously, (|6.2p has a unique solution that satisfies e.g. 
\\Vh{t)\\ < \\VhiO)\\, t € [0,r]. Recah that we have ||u - Pv\\ < Ch'^\\v\\2 for w e iJ^ with v{0) = 0, and 
\\v — Pv\\oc < C'/i^||w||2.cx) for V £ with u(0) — 0. Our aim is to prove an optimal-order error estimate 
for the solution of ()6.2|) . thus extending the result of Dupont, [20], for the periodic initial- value problem for 
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Tjt + Vx — to the case of the initial-boundary- value problem at hand. To do this, we first prove a lemma 
along the lines of the proof of Lemma 5.6. 

Lemma 6.1. Suppose that v is and piecewise on [0, 1] and satisfies v{0) — v'{0) — v"{0) — 0. Then, 
if p = V — Pv and ^ Q is such that 

(C,<^) = (P,0') V0e5L 

we have \\C\\ < C/i2||u(3)||^. 

Proof. We will first show, as in Lemma 5.6, that if — Jj pdx, I < i < N, then Si — 0{h*). Let 

Pv — '^jLi Cj4'j ■ Then, since {Pv,(j)i) = {y,4>i), 1 < i < N, we have Gc — b, where Gy = {(l)j,4>i), 
^ 1^ h i 1^ N , is the N X N tridiagonal matrix with Gu — 2h/3, 1 < i < N — 1, Gnn = h/3, and Gij = h/6 
if \i — j\ = 1, and bi = {v, 1 < i < N. Hence, 



~ J P'^^ ~ J '^'^^ ~ li 1 ^ < i < N , (6.3) 

where 

' ^ if i = 1 
Pvdx ^ { ] ' (6.4) 



^^"'-^+"'^ if 2<i<N. 



The equations of the system Gc — b are 

|(4ci + C2) = 



|(c,_i -|-4c, -hc,+i) = 5i, 2<i<N-l, 

|(cAr_i -I- 2cn) = bN , 

and may be rewritten as 

3 • ^ + ^^^2i+£H) ^ 35^ , 

h{c^-2+c,-i) ^ 4 . hic,^^+cO ^ /^(c,+c.+ o ^ 3^^^_^ ^ ^ 3 < ^ < iV - 1 , 
+ 5 • ^^(^"-^1+^") = 3(6jY_, + 2bN) . 
Therefore, from (|6.4p we conclude that 

r7 = /3, (6.5) 
where 7 = (71, . . . , 7Ar)^, F is the N x N tridiagonal matrix with Fn = 3, F^ = 4, 2 < i < iV — 1, Fatjv — 5 
and T,j 1 if |i - j\ = 1, and /3 = (^1, . . .,I3nV with /3i = 36i, = 3(6.,_i + &»), 2 < i < - 1, and 
/3jv = 3(6jv-i + 26jv)- Hence, from (|6.3p and (I6.5P we have that 

F£ = r, (6.6) 
where e = (ei, . . . ,£Ar)'^, and r = (f'l, . • . tTmY" is given by 

Ti = 3 / wdx -f- / uda; — 3&i , 



Jo Jh 

ri= vdx + 4: / vdx+ / vdx - 3{bi^i + bi) , 2 < i < N - 1 , 

J Xi — 2 "^Xi-i Jxi 

pXM-\ PXM 

rjv = / vdx + ^ I udx — 3(6Ar_i -f 26jv) . 

Jxn~2 Jxn-1 

For ri we have 

ri=3 v{l-f)dx+ v{^-5)dx, 

Jo Jh 

and by our hypothesis on v, using Taylor's theorem with remainder in integral form, we have 

\ri\<Gh^\\v^^^\\L^[oah]- (6.7) 
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As in the proof of Lemma 5.6 we may show, using Taylor's theorem with remainder in integral form, that 



2<i<N . 



(6.8) 



Writing now the matrix F in the form F = 4(/ — E), we see that E is the a N x N tridiagonal matrix with 
PIloo = 1/2. Hence F is invertible, \\T-'^\\oo < 1/2, and (Ill])-(inS give 



Since ((, 



max \e,\ < Ch^v^^^\\oc 

l<i<N 

{p, (p'j), 1 < i < N, we conclude from (|6.9I) that 

max |(C,0,)l<C/i3||«^'^l 

l<'i<N 



(6.9) 



(6.10) 



It is straightforward to check now that Lemma 2.1(i) and (ii) hold for S*^ as well, mutatis mutandis. The 



conclusion of the Lemma follows. 



□ 



Proposition 6.1. Suppose that 7]q E with 77o(0) ~ ?7o(0) — v'oi^) ~ 0- Then, if rjh is the solution of the 
semidiscrete problem (j6'.i?p we have that 

max Wvit) - r,hit)\\ <C{T)h' max hh,oo ■ 

0<t<T 0<t<T 

Proof. Let p ^ T] - Pt], = Pt] - rju. Then, for < t < T 
Putting in this <j) — 9 we obtain, using integration by parts, that 

iTtWmf + ^^'(1,0 = -p{imi,t) + {P.O..). (6.11) 

Now, |p(l,i)| < ||jo||co < Ch?\\'q\\2,oo and our hypotheses on ryo imply that ?/(-,t) e and piecewise for 
< i < T and satisfies dlri(0,t) = for j = 0, 1, 2, < i < T. Therefore, by Lemma 6.1 we have, if 
(C,0i) = (p, 0i)' 1 < * < ^5 that {p,Ox) = {(,9) < C/i^||?7||3,oo||^||- Hence, by the arithmetic-geometric mean 
inequality we obtain from (16. lip that 



i\m)r<ch^r^{t)\\i^^ 



I A 

2 dt 



0<t<T. 



Therefore, Gronwall's Lemma and ()6.2p give the conclusion of the Proposition. □ 

In Table [6T] we show the results of a numerical experiment that we performed to investigate the influence 
of the degree of compatibility of the initial data ryo at a; = on the order of convergence of the L^- error of 
the standard Galerkin semidiscretization of (16. ip using the space for the space discretization and Prjo 





X cxp(a:) 


x"^ cxp(a:) 


x^ cxp(x) 


x^ cxp(a;) 


N 


L'^-error 


order 


L'^-error 


order 


L'^-error 


order 


L'^ -error 


order 


50 


0.9811(-3) 




0.3014(-3) 




0.4786(-3) 




0.6765(-3) 




100 


0.4436(-3) 


1.049 


0.7583(-4) 


1.991 


0.1204(-3) 


1.991 


0.1705(-3) 


1.989 


150 


0.2891(-3) 


1.056 


0.3379(-4) 


1.993 


0.5360(-4) 


1.995 


0.7597(-4) 


1.994 


200 


0.2145(~3) 


1.039 


0.1904(-4) 


1.993 


0.3018(-4) 


1.997 


0.4279(-4) 


1.996 


250 


0.1699(-3) 


1.043 


0.1221(-4) 


1.993 


0.1933(-4) 


1.997 


0.2740(-4) 


1.997 


300 


0.1407(-3) 


1.034 


0.8489(^5) 


1.993 


0.1343(-4) 


1.998 


0.1904(^4) 


1.997 


350 


0.1199(-3) 


1.037 


0.6244(-5) 


1.992 


0.9868(-5) 


1.998 


0.1399(-4) 


1.998 


400 


0.1045(-3) 


1.030 


0.4786(-5) 


1.992 


0.7557(-5) 


1.998 


0.1072(-4) 


1.998 


450 


0.9255(-4) 


1.033 


0.3785(-5) 


1.992 


0.5972(-5) 


1.999 


0.8469(-5) 


1.998 


500 


0.8305(-4) 


1.028 


0.3069(-5) 


1.991 


0.4838(-5) 


1.999 


0.6861(-5) 


1.998 


550 


0.7528(-4) 


1.030 


0.2538(-5) 


1.991 


0.3999(-5) 


1.999 


0.5671(-5) 


1.999 


600 


0.6885(-4) 


1.026 


0.2135(^5) 


1.990 


0.3360(-5) 


1.999 


0.4766(-5) 


1.999 



Table 6.1. errors and convergence rates. Standard Galerkin semidiscretization for 
problem ()6.ip with piecewise linear, continuous functions and four choices of 770(2:). 



as initial data. We considered (|6.ip with T = 0.5 taking successively r]Q{x) equal to a;exp(a;), a;^exp(a;), 
exp(a;) and a;^ exp(x). (The temporal discretization was effected with the Crank- Nicolson implicit scheme 
with time step fc = ft,/3 in all cases; this time step was sufficiently small to ensure that the temporal error 
was about two orders of magnitude less than the spatial error in all cases.) The errors and associated 
rates of convergence of the numerical scheme at t = T = 0.5 for diminishing h = 1 /N are shown in Table [6Tl 
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Recall that if 770(0;) — x'^ exp(x), fc = 1, 2, 3, . . . , the solution r]{x, t) of (|6.ip is in C'^~-^(0, 1) (and piecewise 
C°°) for all t S [0, T]. The LP' orders of convergence turn out to be practically equal to 2 for A: = 4 and A: = 3 
in agreement with the error estimate that was proved. The rate is apparently converging to one in the case 
A: = 1 in which (I6.ip possesses a generalized, piecewise smooth solution, while in the borderline case k — 2 
the convergence rate has not stabilized. (When we increased the number of spatial intervals N, we observed 
that this rate fell steadily and was e.g. equal to 1.9826 at = 1600.) 

6.2. A hyperbolic equation with a variable coefflcient that vanishes at the boundary. Consider 
the problem 

Vt + iui^)x^O, 0<x<l, 0<t<T, 

(6.12) 

ri[x, 0) = ?7o(a^) 1 < a; < 1 , 

where u is a given smooth function that vanishes at x = and at x = 1. No boundary conditions are 
imposed on r]. The problem is well posed and may be solved by the method of characteristics. Uniqueness 
of solutions follows easily: If we multiply the p.d.e. in (|6.12p by 77 and integrate over [0, 1] by parts we see 
that for < t < T 







2dtll'll ' 2 

from which, by Gronwall's Lemma, there follows that 

Um<Choh 0<t<T, 

where C = C(r, maxo<t<T ||u2:||oo)- 

An optimal-order error estimate for the standard Galerkin semidiscretization in S'^ on a uniform mesh 
may be proved for this problem using the methods previously developed in this paper. (This was first pointed 
out in (3j.) The semidiscrete problem consists in finding ry^ : [0,r] — >■ 5^ so that 

(r7;,t,0) + ((u77;,)„0) = O, V0e5L 0<i<T, 

, X (6.13) 
77„(0) = Pr;o, 

where P is the projection onto S^. 

Proposition 6.2. Suppose that u E Cq and that the solution of 1^6.12^ satisfies 77 £ C{0,T; C^). Then 

max ||7/ - r;;i||oo < Ch^ . 

0<t<T 

Proof. Putting 6 = Prj — rjh, p = rj — Prj we have for < < < T that 

(0t,0) + ((u0).,0) + ((up).,0) = O V0e5L (6.14) 
Now, by Lemma 5.1, if for each t £ [0, T], ( £ Sf^ is defined by 

we have that < Ch^. Therefore, putting (p = 6 in (|6.13p and using integration by parts we see that 

||6i|| < Ch^ . 

We conclude that ||r; - 77,,||oo < ||6'||oo + ||p||oo < Ch-^/^\\e\\ + \\p\\oo < Ch\ □ 

This observation may prove useful in the error analysis of first-order problems with advection terms of 
the form {uri)x, where u satisfies Dirichlet boundary conditions. However, as we saw in Section 2 in the 
case of the Boussinesq systems under study the presence of the linear term in the first p.d.e. of (ICB|) or 
(|SCB[) is the obstacle that prevents achieving optimal-order estimates. For example, examining the proof 
of Theorem 2.1 we note that if the term was missing, then the terms {£,x,4') ^-nd {<7xi4') would not be 
present in (|2.21|) and, consequently, by appeal to Lemma 2.2 with w{0) = w{l) = 0, the estimates of the 
terms in (|2.24p would lead to an 0{h^) bound in the right-hand side of (|2.25p . 
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(6.15) 



6.3. Numerical experiments with nonuniform meshes. When we used a quasiuniform mesh in the 
standard Galerkin semidiscretization (|6.2p of the initial-boundary- value problem (j6.1l) . numerical results 
(not shown here) indicating that the error — i]h\\ is of 0{h); of the same order of accuracy is of course 
the upper bound of the error predicted in a straightforward manner by the theory. 

The same differences between the cases of uniform and nonuniform meshes apparently also persist in the 
case of first-order hyperbolic systems. Consider, for example, the following initial-boundary-value problem 
for the nonhomogeneous wave equation written in first-order system form: 

m + ux ^ f , 

0<a;<l, 0<t<T, 

ut + Vx ^ g, 

r/{x, 0) = 'r]o{x) , u{x, 0) = uo{x) , < a; < 1 , 
ri{0,t) = 0, u(l,t)=0, 0<t<T. 

For the purposes of the numerical experiment we took as exact solution of (|6.15p the pair of functions 
r]{x,t) = xeKjp{x{t + 1)), u{x,t) = {x — l)exp(xt) with appropriate initial conditions and right-hand sides. 
We discretized the problem in space by the standard Galerkin method using piecewise linear continuous 
elements on a uniform mesh with meshlength h = and also on a quasiuniform mesh with Ax = 1.6/A^, 
/i2i_i = 0.75Aa;, /i2i = 0.5Ax, 1 < i < N/2. For the temporal discretization we used the classical, fourth- 
order explicit Runge-Kutta scheme with k = h and k = Ax, respectively. As initial values 'r]h{0), Uh(0) we 
took the interpolants of rjo and uq. The errors at T = 0.4 and the corresponding orders of convergence 
are shown, for increasing N, in Table [621[a) in the case of the quasiuniform mesh and in Table [6?2T b) for the 
uniform mesh. The experiment strongly suggests that the errors ore of 0{h?) for a uniform mesh and of 
0{Ax) for quasiuniform. 



N 




order 


-error s{u) 


order 


80 


0.1577(-2) 




0.1607(-2) 




160 


0.7813(-3) 


1.013 


0.8002(-3) 


1.006 


240 


0.5194(-3) 


1.007 


0.5327(-3) 


1.003 


320 


0.3891(-3) 


1.004 


0.3992(-3) 


1.002 


400 


0.3110(-3) 


1.004 


0.3193(-3) 


1.002 


480 


0.2590(-3) 


1.003 


0.2660(-3) 


1.001 


560 


0.2219(-3) 


1.002 


0.2279(-3) 


1.001 


640 


0.1941(-3) 


1.002 


0.1994(-3) 


1.001 


720 


0.1725(-3) 


1.002 


0.1772(-3) 


1.001 



N 


L'^ -error s{ri) 


order 


L'^-errors{u) 


order 


80 


0.8117(-4) 




0.3932(-4) 




160 


0.2029(-4) 


2.000 


0.9825(-5) 


2.001 


240 


0.9021(-5) 


2.000 


0.4369(-5) 


1.999 


320 


0.5074(-5) 


2.000 


0.2457(-5) 


2.000 


400 


0.3248(-5) 


2.000 


0.1573(-5) 


1.999 


480 


0.2255(-5) 


2.000 


0.1092(-5) 


2.000 


560 


0.1657(-5) 


2.000 


0.8026(-6) 


2.000 


640 


0.1269(-5) 


2.000 


0.6145(-6) 


2.000 


720 


0.1002(-5) 


2.000 


0.4856(-6) 


2.000 



(a) 



(b) 



Table 6.2. L^-errors and orders of convergence, problem (|6.15p . Standard Galerkin dis- 
cretization with piecewise linear, continuous functions, (a): quasiuniform mesh (b) : uniform 
mesh. 



As a final remark, we consider the following initial-boundary value problem of a hyperbolic system with 
a viscous term 

0<x< 1, 0<t<T, 

Ut + 'qx- = V> , 

ri{x,0) ^ r]o{x) , u{x , 0) ^ uo{x) , < x < 1 , (6.16) 

u(0,i) = 0, u(l,i) = 0, 0<t<r, 
and discretize it by the standard Galerkin method. We obtain similar results to those of Section 2 in 
which the dispersive term Uxxt is present instead of Uxx- (Cf. Remark 2.3) In a numerical example we 
took a nonhomogeneous version of (I6.16P with exact solution r]{x,t) — exp(2t)(cos(7ra;) -I- a; -I- 2), u{x,t) = 
exp(a;i)(sin(7ra;) -\- x^ — x^) and discretized the problem in space with the standard Galerkin method using 
a uniform mesh with h = 1/N and the quasiuniform mesh of the previous example. We integrated up to 
T = 0.5 with the classical, fourth-order, explicit Runge-Kutta scheme with k = h? /2b in the uniform mesh 
and k = (Aa;)^/25.6 in the nonuniform mesh case. The errors and L'^ convergence rates at T = 0.5 shown in 
Table suggest that the errors for 77 and u are of 0(Ax) and ©(Ax)^, respectively, in the quasiuniform 
mesh case, and of 0{h^^^) and 0{h^) in the uniform case, exactly as in the linear dispersive case. 
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N 










20 


0.2013(-1) 




0.1771(-2) 




40 


0.1004(-1) 


1.003 


0.4434(-3) 


1.998 


60 


0.6696(-2) 


0.999 


0.1971(^3) 


1.999 


80 


0.5023(-2) 


0.999 


0.1109(^3) 


1.999 


100 


0.4019(-2) 


0.999 


0.7099(-4) 


2.000 


120 


0.3350(-2) 


0.999 


0.4930(-4) 


2.000 


140 


0.2872(-2) 


0.999 


0.3622(-4) 


2.000 


160 


0.2513(-2) 


0.999 


0.2774(-4) 


2.000 



(a) 



N 






L ' — GT* OT* S (^U.^ 




20 


0.3605(-2) 




0.1396(-2) 




40 


0.1149(-2) 


1.649 


0.3487(-3) 


2.001 


60 


0.6011(-3) 


1.598 


0.1550(-3) 


2.000 


80 


0.3823(-3) 


1.574 


0.8716(-4) 


2.000 


100 


0.2700(-3) 


1.559 


0.5578(~4) 


2.000 


120 


0.2035(-3) 


1.549 


0.3873(-4) 


2.000 


140 


0.1605(-3) 


1.542 


0.2846(-4) 


2.000 


160 


0.1307(-3) 


1.537 


0.2179(-4) 


2.000 



(b) 



Table 6.3. L^-errors and orders of convergence, problem (I6.16p . Standard Galerkin dis- 
cretization with piecewise linear, continuous functions, (a): quasiuniform mesh (b) : uniform 
mesh. 
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